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ABSTRACT 


8-PSK SIGNALING 
OVER NON-LINEAR 
SATELLITE CHANNELS 
BY 

RUBEN CABALLERO, B.ENG. 


Master of Science in Electrical Engineering 
New Mexico State University 
Las Cruces, New Mexico, 1996 
Dr. Sheila B. Horan, Chair 

Space agencies are under pressure to utilize better bandwidth-efficient communication methods due 
to the actual allocated frequency bands becoming more congested. Also budget reductions is another 
problem that the space agencies must deal with. This budget constraint results in simpler spacecraft 
carrying less communication capabilities and also the reduction in staff to capture data in the earth 
stations. It is then imperative that the most bandwidth efficient communication methods be utilized. 
This thesis presents a study of 8-ary Phase Shift Keying (8PSK) modulation with respect to 
bandwidth, power efficiency, spurious emissions and interference susceptibility over a non-linear 
satellite channel. 
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An end-to-end system performance, including the InterSymbol Interference (ISI) and the Symbol 
Error Rate (SER) as a function of Eb/No on 8PSK modulation was conducted on Signal Processing 
Worksystem (SPW) software installed on a SUN Sparc Station 10 and a Hewlett Packard (HP) Model 
715/100 Unix Station. The simulations were based on the Non-Retum-to-Zero (NRZ) data format. 
The end-to-end system evaluation was performed using idea 1 and non-ideal data with ideal system 
components and three baseband filter types: 5th Order Butterworth, 3rd Order Bessel and Square Root 
Raised Cosine (SRRC), a = 0.25, 0.5 and 1, to observe the effect of pulse shaping on bandwidth and 
SER. 

The simulations show that in-band spurious emissions are greater in number and more evident for 
the Bessel Filter than the Butterworth Filter with a Bandwidth (3dB)-Time product equal to 1 (BT=1). 
With respect to the sideband attenuation, it was found that the values of attenuation for the 3rd Order 
Bessel filter are comparable to the 5th Order Butterworth filter. For SRRC filters with a = 0.25 and 
a = 0.5, the bandwidth is narrower than the Butterworth and Bessel Filters, but the attenuation is less 
at high frequencies. Nonetheless, the absence of spurious emissions is a net advantage. For SRRC 
a = 1, the bandwidth is wider than the two previous SRRC filters and the absence of in-band spurious 
emissions was again noticed. Less sideband attenuation was recorded for this roll-off factor compared 
with a = 0.25 and a = 0.5. For SER, it was found that the Butterworth and Bessel Filters just barely 
meet the threshold of ISI loss < 0.4 dB at SER = 10° . Also the SRRC filters do not meet this 
specification. From these results, it can be stipulated that the 5 th Order Butterworth filter has a small 
advantage compared to the other filters. Overall, it was shown by using baseband filtering that the 
bandwidth utilization can be improved by a factor of 12 to 24 with BT=1 and 8PSK which can 
significantly increase the spectrum utilization. Also a study on the tradeoffs between non-constant 
envelope and bandwidth in a non-linear system. Solid State Power Amplifier (SSPA), is discussed. 
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Chapter 1 


INTRODUCTION AND BACKGROUND 

This thesis gives the results of a study on 8 Level Phase Shift Keying (8PSK) modulation with 
respect to bandwidth, power efficiency, spurious emissions, interference susceptibility and the non- 
constant envelope effect through a non-linear channel. This work was conducted at New Mexico State 
University (NMSU) in the Center for Space Telemetering and Telecommunications Systems in the 
Klipsch School of Electrical and Computer Engineering Department and is supported by a grant from 
the National Aeronautics and Space Administration (NASA) # NAG 5-1491 . 

The first chapter of this report gives a brief summary of the work completed by the Jet Propulsion 
Laboratory (JPL) on various modulation schemes and previous studies performed on non-constant 
envelope signals. In Chapter 2, some theoretical concepts will be explained which are directly related 
to the simulations that were performed for the 8PSK modulation. Finally in Chapter 3, results on the 
simulations for 8PSK will be given. The simulations were performed on a Signal Processing 
Worksystem (SPW) + (software installed on a SUN SPARC 10 Unix Station and Hewlett Packard 
Model 715/100 Unix Station) where the power containment, spurious emissions, symbol error rates 
and non-constant envelope effect on the bandwidth are measured for different types of spectrum 
shaping filters. Conclusions and suggestions for further work are given at the end of this report in 
Chapter 4. 


+ SPW is a registered trademark of COMDISCO systems, a Business Unit of Cadence Design Systems, 
Inc. 919 East Hillsdale Blvd., Foster City, CA 94404. 
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During its 12 th annual meeting (November 1992 in Australia), the Space Frequency Coordination 
Group (SFCG-12) requested the Consultative Committee for Space Data Systems (CCSDS) Radio 
Frequency (RF) and Modulation Subpanel to study and compare various modulation schemes with 
respect to: (a) - bandwidth needed; 

(b) - power efficiency; 

(c) - spurious emissions; and 

(d) - interference susceptibility. 

This study has to be conducted due to the fact that frequency bands are becoming more and more 
congested and space agencies are under constant pressure to reduce costs. To date, four position 
papers have been presented (References [1], [2], [3], and [4]). These reports divided the study into 
three logical phases: 

Phase I (a) - Bandwidth (BW) utilization of various modulation schemes, explores traditional 
modulation schemes that were used in the past and compares them with newer ones that could improve 
the communication channel efficiencies. From this first paper, the following was concluded. The 
traditional modulation methods utilize subcarriers which have some advantages and disadvantages. 
The subcarrier greatly facilitates the separation of data types and also separates the data transmitted 
from the RF carrier. On the other hand, subcarriers significantly increased the spacecraft complexity. 
It also produces additional losses in the modulator and demodulator and occupies a large bandwidth. 
The new modulation techniques (and data formatting) such as PCM/PM/Bi-(j> (Pulse Coded 
Modulated/ Phase Modulation: data are Bi-Phase (Manchester) modulated directly on a residual RF 
carrier), PCM/PM/NR2 (NRZ data are phase modulated directly on a residual RF carrier), BPSK/Bi-tf) 
(data is Bi-Phase (Manchester) modulated on an RF carrier fully suppressing it) and BPSK/NRZ 
(Binary Phase Shift Keying/NRZ data is phase modulated directly on an RF carrier fully suppressing 
it) were examined. 
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It was noted, as expected, that these new types of modulation significantly reduce the amount of 
BW. Therefore use of subcarriers should be limited. A new definition of required bandwidth was 
proposed (required BW is equal to 95% of the corresponding ideally modulated square pulse shaping 
signal) due to the ambiguities of the present definitions. For more information refer to [1]. 

Phase I (b) - A comparison of Quartenary Phase Shift Keying (QPSK), Offset Quartenary Phase 
Shift Keying (OQPSK), BPSK, and Gaussian Minimum Shift Keying (GMSK) [2] is a complimentary 
paper for Phase I (a) which studies modulations that seem more promising for this study (BPSK was 
used as reference). First, it was concluded that by doing some spectral shaping on the data, a 
reduction of the bandwidth can be achieved. The unfiltered data spectrum can be as much as 5 to 10 
times the width of the filtered data spectrum. Also the signal degradation due to power losses, pulse 
distortions and ISI is between 0.2 and 0.4 dB (for QPSK, OQPSK, BPSK, and GMSK) if a matched 
receiver is used. With respect to the pulse shaping filtering, it was mentioned that such filtering at the 
Intermediate Frequency (IF) level would require hardware that would be considered acceptable for RF 
filtering. With respect to the modulation schemes it was noticed that GMSK is a bit better than QPSK 
and filtered BPSK in non-linear channels unless filtering is applied after the power amplifier (PA). In 
the case of pulse shaping filtering, the GMSK and filtered OQPSK offered practically identical 
performance. 

Phase II - Spectrum Shaping [3] is analyzed and simulated. The benefits accruing from the 
spectrum shaping of the transmitted signal are presented. This paper shows that several types of filters 
and location for these filters can be used. Also it was shown that spectrum shaping, with an efficient 
type of modulation, can increase the utilization ratio. This would allow more signals to be placed into 
a given frequency band. The paper also gives the frequency spectra simulations performed on 
PCM/PM/NRZ data on the SPW (COMDISCO) simulation software. Non-ideal and ideal data were 
used with 4 different types of spectrum shaping filters (Butterworth 5th Order, Bessel 3rd Order, 
Raised Cosine (a = 0.25, 0.5 and 1 with NRZ-L data and sampled data), SRRC (a = 1 with NRZ-L 



data and sampled data)). All these simulations were performed using a non-ideal system and at 
passband. Power Spectral plots (Power Containment versus R B ) were generated for these different 
types of filters at the output of the baseband filter, SSPA and 2nd Harmonic filter (refer to [3] for more 
information) and were compared with the case where no baseband filtering is performed. It was 
shown that this baseband filtering in fact improves the frequency band Utilization Ratio. Studies were 
also performed to determine the filter amplitude response and ISI of these filters. The Raised Cosine 
filter was found to be amplitude unstable at high data transition rates therefore it was not included in 
Phase III. 

As mentioned in the paper of Phase II, Phase III should apply the techniques (Power Containment 
Plots, Utilization Ratio...) described in Phase II to each of the modulation schemes identified in Phase 

I, i-e., 

(a) - PCM/PM/NRZ (NRZ Data Modulated directly on RF Carrier); 

(b) - PCM/PM/Bi-<j) (Split-Phase Data Modulated directly on RF carrier) ; 

(c) - BPSK/NRZ (NRZ Data modulated on suppressed carrier); 

(d) - BPSK/ Bi-(f) (Split-Phase Data Modulated on suppressed RF Carrier); 

(e) - QPSK (NRZ Data Quadri-Phase Modulated RF Carrier); 

(f) - OQPSK (NRZ Data Offset QPSK Modulated on RF Carrier); 

(g) - 8PSK (NRZ Data 8-Phase Shift Keyed Modulated on RF Carrier); 

(h) - MSK (NRZ Data Minimum Shift Keyed on RF Carrier); and 

(I) - GMSK (NRZ Data Gaussian Minimum Shift Keyed on RF Carrier). 

End-to-end system performance, including the ISI and the SER as function of Eb/No should be 
determined. Therefore at the end of Phase III, the modulation techniques providing the most efficient 
frequency spectrum utilization should become clear. 
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Phase III (a) - End-to-End System Performance for PCM/PM/NRZ. Reference [4] contains the 
simulations performed on SPW software. The simulations were based on the PCM/PM/NRZ 
modulation format. The end-to-end performance was performed using non-ideal data, non-ideal 
system components and three baseband filter types 

(a) - Butterworth 5 th Order; 

(b) - Bessel 3rd Order; and 

(c) -SRRC (a= 1). 

Again note that the Raised Cosine filters were not included in Phase III since they were discovered 
to be amplitude unstable at high data transition rates. The idea is to have a design utilizing the 
narrowest possible transmitted RF Spectrum and yet be consistent with reasonable system losses, (i.e. 
maximize bandwidth efficiency with acceptable ISI loss). The optimum filter is found by selecting the 
filter that provides the lowest Bandwidth-Time product (BT) where T is the duration of one symbol 
(for PCM/PM/NRZ the duration of one bit is equal to the duration of one symbol or 1/R B where R B is 
the bit rate) . From these simulations the following conclusions were derived for PCM/PM/NRZ 
modulation. The simulations show that a practical system can be built with filtering ISI losses below 
0.4 dB and the total end-to-end performance system losses below 0.75 dB at SER = 10°. Also 
excellent utilization ratios are achieved based on -50 dB sideband attenuation. The following table 
summarizes the results for the different types of baseband spectrum filters used: 


Table 1.1 - Summary of PCM/PM/NRZ Modulation Simulation 


PARAMETER 

BUTTERWORTH 

BESSEL 

SRRC 

BT or a 

2 

2 

1 

ISI LOSS (dB) 

0.28 

0.37 

0.19 

TOTAL SYSTEM LOSS (dB) 

0.64 

0.73 

0.55 

UTILIZATION RATIO, p 

8.5 

7.4 

8.1 
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As shown in this table, significant bandwidth savings (Utilization Ratio) is possible by selecting 
the proper baseband filter. Different modulation types will also affect the bandwidth savings. 
Therefore, future work will involve the same measurements as described above for the different types 
of modulation. It will then be possible to identify the most efficient modulation type(s) based upon the 
application, the system performance and the preferred frequency band. In this case , this report 
contains the simulations performed on an 8PSK signal being filtered (spectrum shaping) by the three 
filters mentioned earlier. 

To complete this work, the relationship between the non-constant envelope introduced by pulse 
shaping the 8PSK signal and the bandwidth of the filter was to be investigated. Therefore the 
following question had to be answered: 

How much variation is introduced by the non-constant envelope 
as the bandwidth of the filter is made smaller? 

A literature search was conducted to see if any kind of investigation had been performed with 
respect to an 8PSK non-constant envelope signal going through a non-linear channel (amplifier). The 
analysis of the non-linear distortion caused by the amplifier, more specifically for a non-constant 
envelope, is very complicated. Although many analyses have been published in the literature, none 
has been entirely applicable here. The research conducted in [5] shows how a bandlimited spectrum 
filter and the shape of the frequency pulse affect the error probability of a frequency shift keying 
(FSK) signal with differential phase detection. The investigation was conducted using a Butterworth 
filter (order = 4 and 3 for the transmitter and receiver respectively) and the frequency shaping pulse is 
raised cosine or rectangular. This investigation was then performed for a FSK signal as opposed to a 
PSK signal. In [6], an experiment more related to the one described in this thesis was performed. 
This paper investigates the effect of filtering PSK, FSK and MSK signals in a theoretical way. A 
Butterworth Filter (4th Order) is used and the ISI and BER are measured. 
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However, this paper does not consider an amplifier in their study. In [7] , the performance of a 
satellite channel using a software package (developed at Telecom Research Laboratories) was 
measured. 

The results of this experiment (using the software) were then validated by experiments on a 
hardware simulator. It was shown that both results from the software and actual measurements are 
very similar (minor differences due to the ideal demodulator used in the software package). Even 
though the satellite channel in this paper uses spectrum shaping filters and non-linearities (TWT 
Amplifier), there was no experiment conducted to measure the effect of non-constant envelope 
through this non-linear system. 

Although other similar studies were performed on pulse shaping, non-constant envelope and non- 
linear channels, none of them involved the measurement of the effect of a non-constant envelope 
signal (8PSK for this case) through a non-linear channel (SSPA). The following thesis gives the 
results on simulations performed on an 8PSK bandlimited signal ( bandlimited with the three filters 
enumerated earlier) going through a non-linear amplifier. The effect of the non-constant envelope 
through this channel will be demonstrated. 
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Chapter 2 


THEORY 


This section will briefly describe some of the concepts that were used for the simulations performed 
on the SPW simulation software (from COMDISCO) that is installed on a SUN Sparc Station 10 and 
HP Model 715/100 Unix Station in the Center for Space Telemetering and Telecommunications 
Systems in the Klipsch School of Electrical and Computer Engineering Department at NMSU. 

2. 1 Eight (8) Level Phase Shift Keying (8PSK) 

8 level PSK is a multi-phase signaling technique meaning that 8 modulation phases (with constant 
amplitude) are produced at the transmitter output. Figure 2.1 shows how an 8-phase signal can be 
generated from a serial binary input stream using a 3 bit Serial-to-Parallel Converter (SPC). In this 
case, since the data rate is R B bits/sec, then the Symbol Rate (baud), R s , of the multi-phase signal will 
be 

Rs = ReA where X =log 2 8 = 3 thus R s = R B /3 (symbols/sec). 


Binary Input 
Data Rate = R 


^ 

3 bits 

3 bits in parallel 
. — — 

TRANSMITTER 


B bits/sec 

SPC 

Rs = Rg/3 (symbols/sec) 

1 


OUTPUT 

8PSK 


Figure 2.1 - 8PSK Digital Transmission System 


A plot of the complex envelope of the 8PSK signal, g(t) = A c e je(t) , is shown in Figure 2.2. This 
figure is usually referred as the 8PSK constellation. 
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IMAGINARY 

(QUADRATURE) 

A 

Oil 001 


010 • 


110 • 


y. 


000 

Ag, : 

X, 


• 100 


REAL 

(IN PHASE) 


• • 

111 101 


Figure 2.2 - 8PSK Signal Constellation 

As shown in Figure 2.2, each point on the constellation corresponds to one of the 8 multi-phase 
values (or symbols) which corresponds to one of the phases that 0 is allowed to have. In the above 
figure each point corresponds to a symbol where the PSK phases are 

0, =22.5° for symbol 000 (decision region: 0° to 45°) 

0 2 = 67.5° for symbol 001 (decision region: 45° to 90°) 

0 3 = 1 12.5° for symbol 01 1 (decision region: 90° to 135°) 

0 4 = 157.5° for symbol 010 (decision region: 135° to 180°) 

0 5 = 202.5° for symbol 1 10 (decision region: 180° to 225°) 

0 6 = 247.5° for symbol 1 1 1 (decision region: 225° to 270°) 

0 7 = 292.5° for symbol 101 (decision region: 270° to 315°) 

0 g = 337.5° for symbol 100 (decision region: 315° to 360°) 
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Another way of expressing 8PSK is by using two orthogonal carriers modulated by x and y 
components of the complex envelope (no phase modulator is used in this case) thus 

g(t) = A c e ,ew = x(t) + jy(t) 

where x, = A c cos 6, and yi = A c sin 0; where i = 1 to 8 for 8PSK as shown in Figure 2.2. 

For the bandpass physical signal, the waveform can then be represented by 

v(t) = Re{g(t)e’ <0c '} where co c = 2jtf c , f c is the carrier frequency in hertz (Hz). 

For the simulations performed on SPW, the 8PSK modulator in the library was slightly modified to 
be able to implement the non-ideal data. By changing some of the constants, this 8PSK modulator can 
become a BPSK or QPSK modulator. The complex envelope form of 8-PSK was used to develop this 
baseband modulator. Figure 2.3 shows the block diagram used in SPW. The “Control Block” in the 
figure connects one of its two inputs (xl or x2) to the output depending on the state of the logical 
“control” input. When the control input is 0, xl is used; when the control input is 1, x2 is used. 



Figure 2.3 - Block Diagram for the BPSK, QPSK and 8PSK Modulator on SPW 
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The values of Q are given in Table 2.1 and vary depending on the desired modulator (BPSK, 
QPSK or 8-PSK) 

Table 2.1 - Values of Q for the BPSK, QPSK, 8PSK Simulation Model 


MODULATION 



VALUES 

OF CONSTANT “C” 



TYPE 

Cl 

C2 

C3 

C4 

C5 

C6 

C7 

BPSK 

0.0 

0.0 

0.0 

2n/2 

0.0 

0.0 

42 

QPSK 

2.0 

0.0 

0.5 

2tt/4 

0.0 

0.0 

42 

8PSK 

2.0 

4.0 

0.5 

2ti/8 

0.0 

0.0 

42 


2.2 NRZ-L Ideal And Non-Ideal Data 

The binary input or data format that is used throughout these simulations is referred as NRZ-L 
which stands for Non-Retum-to-Zero Logic level. NRZ-L has the logic 0 or 1 represented by the 
appropriate voltage level for the duration of the bit period. The following convention will be used: a 
logic 0 will be represented by a -A Volt (V) signal level and a logic 1 will be represented by a logic 
+A V signal level (which is the convention used in these simulations). Figure 2.4 shows this digital 
data format. 

BINARY DATA 

NRZ-L 1 10 1 0 0 ; 1 



Bit Period * T B 

Figure 2.4 - NRZ-L Digital Data Format 
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All the simulations were done using ideal and non-ideal data. Ideal data are defined as having a 
perfect symmetry, i.e., the duration of a digit one (1) is equal to the time duration of a digit zero (0) 
and it also has a perfect data balance, i.e., the probability of getting a zero is equal to the probability of 
getting a 1 (Pr(0) = Pr(l) = 0.5 or 50%). For non-ideal data these two conditions, data symmetry and 
data balance, are not respected. The CCSDS limits are ± 2% for data asymmetry and a data imbalance 
of 0.45 (probability of getting a 1 vs probability of a 0) as mentioned in [3], The data asymmetry can 
be increased by the stray capacitance in spacecraft wiring and data imbalance can be produced by long 
runs of Is and 0s in the random data. If non-ideal data are present (i.e., the mean value or expected 
value, m a , of the signal is not equal to 0), the Power Spectrum Density (PSD) of the digital signal 
will then consist not only of a continuous spectrum that depends on the pulse-shape spectrum of the 
signal data (rectangular pulse for NRZ-L), but will also contain spectral lines (delta functions) spaced 
at approximately the harmonics of the symbol rate , R s . The PSD equation of a baseband digital 
signal for the case of uncorrelated data, derived in [8] shows these results: 


PSD x (f)= C , 2 ./V|F(/f 


Hm. R s ) 2 - j|F(nR s )|5(/-nR s ) 

n=- oc 


Continuous spectrum 


Discrete Spectrum 


where m a is the mean of the data 

a a 2 is the variance of the data 
R s = 1/T S = Baud or Symbol Rate 

F(f) is the pulse shape spectrum (sine 2 shape for NRZ-L data). 
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2.3 Intersymbol Interference (ISI) 

ISI can be introduced into a system by improperly filtering pulses as they pass through a 
communication system. As it is known in communications, the absolute BW of flat-top multilevel 
pulses is infinity, therefore since the channel is always bandlimited, the pulse waveforms passing 
through it will be dispersed or spread. This spreading will exceed the time duration of the pulse or 
symbol and will overlap with the next symbol. This interference is called ISI which will cause system 
degradation or higher error rates. Figure 2.5 shows two examples of ISI. 


Signal sent 
in the channel 

1 0 0 0 


Impulse Response of the signal 
through the channel 



Signal sent 
in the channel 

1 0 ! 1 1 



Time (t) 


Impulse Response of the signal 
through the channel 


ISI 



0 sec 


Bit Period = T B 


Bit Period = T n 

Sampling 
Times 



Figure 2.5 - Intersymbol Interference (ISI) 
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In 1928, Nyquist discovered three different methods to eliminate ISI by having rounded tops 
pulses instead of flat tops [8]. The BW can then be restricted without introducing ISI. One of these 
methods includes the utilization of Raised Cosine - Rolloff Filtering that will be described below. 


2.4 Spectrum Shaping And Filters 

Spectrum Shaping is a technique that allows the signal energy to be concentrated near the carrier 
and therefore reduce the power contained in the sidelobes (baseband and RF). This is done by filtering 
the signal. Such a technique can also yield a significant savings in the required BW. Nonetheless, 
this type of filtering causes loss of energy (if filtering is done at the RF level) and ISI or pulse 
distortion (if the BW limitation is too restrictive since the pulses are not rectangular any more). These 
losses can be reduced if the receiver compensates for the distortions (i.e., if a matched filter is used). 
Thus baseband filtering is effective in limiting spectrum width. 

The filters that were used for the simulations are Butterworth, Bessel and Raised Cosine. They will 
be described in the following paragraphs. 

2.4.1 Butterworth Filters 

The Butterworth family of filters is known to have a “maximally” flat response in the passband 
region [9]. In other words the amplitude response in the passband region comes very close to a 
perfectly flat response (at a given order). This is why Butterworth Filters are commonly used. Figure 
2.6 shows the amplitude response for different orders of Butterworth filters. Note again how flat the 
amplitude is in the passband region. All orders of this filter pass through the 3 dB cut-off point. Also 
as the order of the filter increases, the sharpness of the slope after the 3dB cutoff point also increases. 
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Butterworth Filter fc(3dB) = 1 Hz& order=2, 3 ,4,5,67.8 



Figure 2.6 - Amplitude Response of Various Orders of Butterworth Filters 

2.4.2 Bessel Filters 

The main characteristic of the Bessel family of filters is that they have an approximately linear 
phase in the passband region. This linear phase will have a tendency to prevent dispersion of the 
signal and therefore can be good for digital pulses. Figure 2.7 shows the amplitude response of 
various orders of Bessel Filters. Note that the amplitude is not as flat in the passband region as for the 
Butterworth filter. Also the transition from the passband to the stopband region is not as rapid as for 
the Butterworth Filters for the various orders of filters. 
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Bessel Frrtrer fc(3dB) = 1 Hz & n=2, 3,4,5,67,8 



Frequency (Hz) 


Figure 2.7 - Amplitude Response of Various Orders of Bessel Filters 
For the simulations performed in this paper, a 5th Order Butterworth Filter and a 3rd Order Bessel 
Filter were selected. These orders of filters are the same as those selected for comparison during Phase 
III of this study. 

2.4.3 Raised Cosine Filters 

As mentioned earlier, if the Raised Cosine filters are used the ISI will be eliminated in a linear 
channel. The idea behind the Raised Cosine Filter is that the time signal goes through zero at adjacent 
sampling points therefore eliminating the interference of other symbols. Also it has the important 
property of having an envelope that decays fast enough so that clock phase jitter at the sampling times 
do not introduce to much 1ST 
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The transfer function of the Raised Cosine Filter is 


) 

0+g) 

2T 


where the zeros will occur at t=nT (T is the sampling interval) and a is the rolloff factor or the excess 
bandwidth over the Nyquist Band Filter which can be varied from 0 to 1. The rolloff factor, a, 
determines the bandwidth of a Raised Cosine Filter. For a = 0, the filter’s transfer function 
approximates that of a “Brick Wall” filter with bandwidth 1/(2T) while an a=l yields a sinusoidal 
transfer function having a total width of 1/T. Figure 2.8 shows the frequency response of the Raised 
Cosine Function with different rolloff factors (a = 0 , 0.25, 0.5 and 1) and T =1 seconds. 


Transfer Function of Raised Cosine Filter with r = 0, 0.25, 0.5 and 1 



X(J) = 


\T, 

L 

2 v 

[0, 


1 - sin 


(2n\f\T-%) 


T\ 


0 <|/|< 

(1-a) 


2a 


2T 


(i-a 

2T 

£|/|* 


l/i 


(i+g) 

2T 


Figure 2.8 - Frequency Response: Raised Cosine Characteristics with Different a 
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Figure 2.9 gives the time domain response for different rolloff factors. 


Time Response of Raised Cosine Filter with r = 0, 0.25, 0.5 and 1 



Time (t) 

Figure 2.9 - Time Response of Raised Cosine Filters for Different a 
From Figure 2.9, it can be seen that if the symbols are sampled at the nulls of the sine function 
then the symbol interference from the previous and next symbols will not exist. This is why the 
Raised Cosine filters are used since it is a class of filter that satisfies Nyquisf s first criterion (see [8]). 
To produce a zero IS1 channel the impulse response shown above has to be created. This can be done 
by utilizing a Square Root Raised Cosine filter in both the transmitter and the receiver. The overall 
filtering characteristics of the channel will be a Nyquist Channel or a Raised Cosine Filter: 


X(f) = T(f) • H(f) • R(f) = < 

Y>0<|f|<T 


0, elsewhere 
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where T(f) is the frequency response of the transmitter; 

H(f) is the frequency response of the channel; 

R(f) is the frequency response of the receiver; and 
X(f) is the overall frequency response of the system. 

If the channel has a large bandwidth compared to the transmitter and receiver frequency response 
then 

W) = r(/)j?(/) = RC(f) 

where RC(f) is the frequency response of the Raised Cosine Filter. Thus 


TV) = RV) = JRCV) = SRRC(f) 


where SRRC(f) is a Square Root Raised Cosine Filter for both the transmitter and the 
receiver which maximizes the data transmission efficiency (both the transmitter and receiver are 
matched). Therefore, the channel will be a Nyquist Channel (or Raised Cosine Function) with zero ISI 
for optimum data detection. This does not apply for a non-linear channel as will be shown when a 
non-linear device such as a Solid State Power Amplifier (SSPA) is present in the channel. 

As stated in [3]: “Raised Cosine Filters were selected for evaluation because the linearity of their 
phase-frequency relationship should help to eliminate the ringing found in the Butterworth filters at the 
cut-off frequency. Their comparatively narrow bandwidth, combined with a smooth response, should 
provide a signal which concentrates most of the data sideband energy near the main lobe significantly 
attenuating the sidebands. Such filters are commonly employed to pack a multiplicity of signals in a 
confined frequency band.” 
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2.5 Solid State Power Amplifier (SSPA ) 


Amplifiers are classified into two main categories: linear and non-linear circuits. Devices such as 
Traveling-Wave Tubes (TWT) and Solid State Power Amplifier’s (SSPA) can be linear amplifiers 
when operated well below the saturation level. If this drive level is increased, the efficiency, Le., RF 
output/Direct Current (DC) input is improved, but the amplifier then becomes a non-linear device. 
Therefore constant envelope signals such as PSK are preferred so that the distortion due to the non- 
linear amplifier will be minimized. In this case, the power amplifier will be used in its non-linear 
region to maximize the signal power. Figure 2.10 shows the SSPA, European Space Agency (ESA) 10 
Watts, Magnitude Characteristics Curve. 


SSPA 10 Watts Amplifier 



Input Level Power (dB) 


Figure 2.10 - SSPA, ESA 10 Watts, Magnitude Characteristic Curve 
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The non-linearities that occur in the nonlinear region of this amplifier will cause non-linear 
distortion and Amplitude Modulation-to-Phase Modulation (AM-to-PM) and AM-to-AM conversion 
effects (linear channels are channels without AM-AM and AM-PM conversions). As mentioned in 
[8], the analysis of these non-linearities are very complicated. 

In this study the non-linear channel will be considered (the power amplifier will be operated in its 
saturation region). Even if SRRC filters are used for the transmitter and receiver, ISI-ffee sample 
points no longer exist because of the distortion introduced by the system. This will be shown and 
discussed in Chapter 3 - 8PSK Simulations. Thus, the objective of this study is not only to find a 
bandwidth efficient communications system which can increase frequency band-utilization, but also, to 
identify an implementation that can be realized. 

NMSU used the SSPA model for their simulation which is based upon specifications provided by 
the European Space Agency (ESA) for their 10 Watts, solid state, S-band power amplifier. This power 
amplifier was selected since all the simulations performed in Phases II and III (refer to [3], [4]) were 
done using this amplifier. 

2.6 Additive White Gaussian Noise (AWGN) 

In the channel being simulated, white gaussian noise was added to the signal. In AWGN, the term 
additive means that the noise is added to the signal. Also by Gaussian it is meant that the random 
process representing the noise, n(t), is Gaussian and has a Power Spectral Density (PSD) of 


PSD n (f)= 


No/2, \f\ < B 
0 , otherwise 


where B is a positive constant which represents the BW and the noise is bandlimited as long as B is 
finite. If B-» oo, the noise is completely white, i.e., all the frequencies are present. 
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Gaussian processes g,ve a good approbation of naturally occurring behaviors. Addition of the 
noise, n(t), will give r(t) = s(t) + n(t) where 

- n(t) is the AWGN with variance = a 2 = No/2 and mean = m n = 0; 

- s(t) is the signal; and 

- r(t) is the resulting signal with the noise added. 


This resulting signal (random process) will have a Gaussian distribution with mean - s(« whete c, 
is the sampling tune and a variance o,= - No/2, figure 2. 1 1, shows a block diagram of how the notse 
was added to the initial signal on SPW. Figure A.4 in Appendix A contains a more derailed diagram 
of the circuit used in SPW. As shown in Figure 2.1,, the variance ( and therefore the standard 
deviation ) is calculated by using the average complex power of the signal and Es/No (dB) whtch was 
entered by the user. The output of this block gives the input signal added whh a Gaussian Noise. 


Controlled Es/N 0 
(not in dB) 



Signal 

plus 

Noise 


Figure 2.11- Block Diagram for the AWGN 
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In our simulations, a variation of the integrate-and-dump filter was used. This was possible since 
the digital signaling format has a rectangular bit shape. For the integrate-and-dump filter, the digital 
input signal plus the noise is integrated over one symbol period T and the output of the integrator is 
dumped at the end of the symbol period which is the value that is required. Therefore for proper 
operation of this optimum filter, an external clocking signal called bit sync is required. 
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Chapter 3 


8PSK SIMULATIONS 


3.1 Simulation Procedures 


To simulate the non linear satellite link, a Signal Processing Worksystem (SPW) which is a 
COMDISCO simulation software was utilized on a SUN Sparc Station 10 and HP Model 715/100 


Unix Station. The following diagram shows the different blocks that were simulated on the 
workstation. 



Figure 3.1 - Simulation Block Diagram 
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First a Data Source or Modulator was used to produce the 8PSK signal. The data source also 
contained a block that could produce data asymmetry on the data. All simulations were produced at 
baseband using the complex envelope signal representation. The results do not vary if the simulations 
are done in baseband or passband. The advantage of going to baseband is the computer execution 
time, i.e., it takes longer to simulate at passband since the sampling frequency must be at least twice 
the Nyquist rate. The next block after the modulator is the baseband spectrum shaping filter. The three 
types of filters that were used for these simulations were 5th Order Butterworth, 3rd Order Bessel, and 
finally Square Root Raised Cosine (SRRC) with Sampled Data and Roll-Off factors of 0.25, 0.5 and 1 . 

The next blocks include the SSPA, a bandlimiting filter and the noise (channel). Simulations were 
performed using the SSPA at the saturation level to maximize the power. Simulations at saturation 
were then followed by simulations with the SSPA backed off by -10 dB. This amplifier was followed 
by a 2nd Harmonic filter (4th Order Butterworth with a bandwidth of ± 20Rs) which reduces the 
interference between different channels. This bandlimiting filter is followed by the variable AWGN 
(variable to be able to control the Es/No values). 

These blocks were followed by the Receiver, Symbol Synchronizer and finally the Error Estimator. 
A matched filter (sliding integrator) was used as the receiver for the Butterworth and Bessel Filters. 
The matched filter was matched to the NRZ-L baseband data. For the SRRC, the receiver was also a 
SRRC to minimize the ISI. The SRRC will minimize the IS1 and the matched filter is an optimum 
filter which optimizes the Signal-to-Noise ratio (SNR). The synchronizer that was used consisted in a 
Delay and Phase Meter that would correlate, as best as it could, the initial data with the data that went 
through the channel. 

After these two signals went through the phase and delay meter, the initial data were delayed to be 
synchronized with the distorted data (the delay between the initial and distorted data were caused by 
the filters). After the data were synchronized (the best correlation was found), an Error Rate 
Estimator was used to measure the differences in symbols and calculate the Symbol Error Rate (SER). 
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When the Power Containments were needed a Spectrum Analyzer from SPW was placed after the 
2nd Harmonic Filter. 

The following variables and systems were implemented in SPW for the simulations: 

(a) - SPW Simulation Variables: 

i - Bit rate - 1 bps (R s = symbol rate = ( 1/3) symbols/s for 8PSK) 

ii - Sample Rate = f s = 16 samples/sec 

iii - Data: NRZ-L 

iv - Carrier Frequency = 0Hz (Baseband Simulations) 

(b) - Transmitting System: 

i - Data Generator: Ideal and Non-Ideal Data: 

- data asymmetry = 2% 

- data imbalance = 0.45 

ii - Baseband filters (do not include resistive and reactive losses): 

- Butterworth 5th Order (BT =1,2,3); 

- Bessel 3rd Order (BT= 1,2,3) 

- SRRC (a =0.25, 0.5 and 1) with 256 taps 

Note - BT represents the product of the bandwidth (3 dB single-sided) 
with the symbol time; thus if BT = 1 it means that 
B = (1/3) since T = (1/Rs) = (l/(l/3)) = 3. 

- Also 256 taps where chosen because we will have: 

(256 taps/ 16 samples/sec)/2 = 8 symbols interfering with 
the one being sampled. It was noted that using more 
taps on the Infinite Impulse Response (IIR) filter did not 
affect the output. 
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iii - Power Amplifier based on: 

- European Space Agency (ESA) 10- Watt 
Solid State Power Amplifier (SSPA); 

- 4th Order Butterworth (cutoff at ± 20 R s ) - 2nd Harmonic Filter 

- protects user’s of other bands. 

(c) - Receiving System: 

i - Matched Filter: Sliding Integrator 

ii - Delay and Phase Meter (synchronization); and 

iii - Error Rate Estimator. 

The following simulations were performed with the system described above 

(a) - Power Containment and Spurious Emissions; 

(b) - End-to-End System Performance : Symbol Error Rate (SER); and 

(c) - Non-Constant Envelope Simulations (a variation of the circuit shown in Figure 3.1 was 

utilized). 

3.2 Power Containment and Spurious Emissions Simulations 

Table 3.1, summarizes the simulations performed for 8PSK ( the numbers in parenthesis, (1) etc. 


indicate individual runs that were performed): 



Table 3.1 - Power Spectrum Simulations 


PLOTS 

LOCATION 

( 1 ) - Modulator Output 

(2) - Output of the SSPA 

(3) - 2 nd Harmonic Filter Output 

(1) - Modulator Output 

(2) -Output of the SSPA 

(3) - 2 nd Harmonic Filter Outnut 

(1) - Spectrum Filter Output 

(2) -Output of the SSPA 

(3) - 2 nd Harmonic Filter OutDut 

(1) - Spectrum Filter Output 

(2) - Output of the SSPA 

(3) - 2 nd Harmonic Filter Outnut 

(1) - Spectrum Filter Output 

(2) -Output of the SSPA 

(3) - 2 nd Harmonic Filter Outnut 

( 1 ) - Spectrum Filter Output 

(2) -Output of the SSPA 

(3) - 2 nd Harmonic Filter OutDut 

( 1 ) - Spectrum Filter Output 

(2) - Output of the SSPA 

(3) - 2 nd Harmonic Filter Outnut 

(1) - Spectrum Filter Output 

(2) - Output of the SSPA 

(3) - 2 nd Harmonic Filter OutDut 

(1) - Spectrum Filter Output 

(2) -Output of the SSPA 
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(1) - Spectrum Filter Output 
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(3) - 2 nd Harmonic Filter Output 
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SIMULATION fl) - UNFILTERED BASEBAND IDEAL DATA 


Figure 3.2a shows the power spectrum of the NRZ-L Ideal data using 8PSK. The abscissa is the 
frequency with respect to the bit rate R B thus the first lobe for 8PSK ends at ± R B /3 (R B /3 represents 
the baud rate or symbol rate R s ) which was to be expected. The following table shows the comparison 
of the theoretical and simulated amplitudes for 8PSK using SPW with Ideal data. The values in the 
simulated amplitude column were taken from Figure 3.2a by moving the cursor on the peak of each 
lobe at the frequencies indicated; the accuracy of the readings of the simulated amplitude was 
approximated to be around ± 1 .5 dB. It is to be noted that the results of the simulations are very close 
to theory. The following formula was used to calculate the theoretical power spectrum [8]: 


PSD(f) = 


J sin(izfil b )\ 




where M-2* is the number of points in the signal constellation ( in this case X= 3 for 8PSK) and 
the bit rate is R B = 1/T B (1 bit/sec in this case). 


Table 3.2 - Comparison of Theoretical and Simulated Amplitudes for 8PSK using SPW (Ideal Data) 


OFFSET FROM 
CENTER FREQUENCY 
(±Rb) 

THEORETICAL 

AMPLITUDE 

(dB) 

SIMULATED 

AMPLITUDE 

(dB) 

0 

0 

0 

0.5 

-13.46 

-13.5 | 

0.8333 

-17.90 

’ -18.8 

1.16667 

-20.82 

-21.3 

1.5 

r -23.01 

-23.9 

1.8333 

-24.75 

r~ -25.6 

2.1667 

-26.20 

-27.1 

2.5 

-27.44 

-28.3 

2.8333 

-28.53 

-29.2 

3.1667 

-29.50 

-30.9 

3.5 

-30.37 

-31.0 

3.8333 

-31.16 

-31.9 

4.1667 

-31.88 

-32.5 

4.5 1 

-32.55 

-33.2 1 

4.8333 

-33.17 

-33.7 

5.1667 

-33.75 

-34.5 [ 
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Figure 3.2 - Power Spectra of Unfiltered 8PSK Baseband Ideal Data 
using SSPA at Saturation Level 




These data and Figure 3.2a will be used as reference for the next figures. To produce Figure 3.2a, a 
sampling frequency, f s , of 250 samples/sec was chosen. The sampling frequency was chosen much 
higher than the Nyquist rate (2 times the highest frequency) to be able to have a spectrum going from 
-125 R b to + 125R B . Also 300 000 samples (to have 1200 symbols, i.e., 1 symbol/sec *3 00 000 
samples/ 250 samples/sec) were taken therefore an FFT of N=262144 (power of 2) was utilized. 
Hence the frequency resolution, f , , [10] of the plots is given by: 

frequency resolution = f, = ^ =9.537x10"“ Hz 

Figure 3.2b shows the same results as Figure 3.2a but the frequency abscissa was set to ± 50 R B . 
Figure 3.2c shows the output spectrum of the SSPA with ideal data. The SSPA was used at its 
saturation level (0 dB). Simulations were later done with the SSPA set to 10 dB away from saturation 
(10 dB Backoff) The output of the SSPA is identical to it’s input due to the constant envelope of the 
input signal. Figure 3.2d is the same as Figure 3.2c but at ± 50 Rg instead of at ±5 R B . Finally 
Figure 3.2e shows the output of the 4th Order Butterworth (± 20 Re) 2nd Harmonic filter which has no 
effect in the region of ± 5 R B . This filter is used to reduce unwanted emissions on other bands. Figure 
->.2f is the same as Figure 3.2e but ± 150 R B was used in the abscissa to see the effect on the other 
bands. 
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SIMULATION (2) - UNFILTERED BASEBAND NON-IDEAL DATA 


For this simulation, the data symmetry was set to 2% and the probability of zero = 0.45. Figure 
3.3a shows the unfiltered baseband data (Non-Ideal Data) at the output of the 8PSK modulator. The 
presence of spurs at approximately f = n*R s (where n=0,l,2..) are due to the probability of zero which 
is not equiprobable, i.e., equal to 0.5, and to the data asymmetry. The spurs are more evident in Figure 

3.3b. 

Figure 3.3c and 3.3d show the output of the SSPA (at saturation level) and as mentioned earlier 
there is no effect or difference with the input data since it is a constant envelope signal (no spectrum 
shaping filter present) going through the amplifier. Figure 3.3e and 3.3f show the output of the 4th 
Order Butterworth filter at ± 5R B and ± 1 50R B . 

Simulations (1) and (2) will be used as references for the spectrum shaping filter simulations. 
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SIMULATION (3) AND (4) - 5TH ORDER BUTTER WORTH FILTERED BASEBAND DATA 

(IDEAL AND NON-IDEAL DATA) 

For this simulation, the sampling frequency (f s ) was set to 250 samples/sec, the symbol rate (R s ) = 
(1/3) symbols/sec, the number of samples set to 300 000 and the number of FFT points (N) = 524 288. 
The BT product was set to 1 thus the bandwidth B = R s and T = 1/R S which is the symbol interval. 
For the non-ideal data, the data asymmetry was again set to 2% and the probability of zero was set to 
0.45. 

Figure 3.4a and 3.5a show the result of the baseband data going through the 5th Order Butterworth 
Filter for ideal and non-ideal data, respectively. Note that at approximately ± 2R S the filter attenuates 
the data sidebands by approximately 40 dB. Figures 3.4c and 3.5c show the output of the SSPA for 
the ideal and non-ideal data. Note that due to the non-constant envelope of the data (since some 
spectrum shaping was performed), the output of the SSPA is quite different from Figure 3.2c and 
3.3c. In fact, it seems like the SSPA is trying to recreate the sidelobes that were attenuated by the 
filter. Also, the spurious emissions encountered in Figure 3.3d for the non-ideal case are not present or 
are fairly attenuated. In [3] for PCM/PM/NRZ it was reported that the Butterworth Baseband Filter 
introduced large in-band spurious emissions which does not seem to appear here. Comparing Figures 
3.4d and 3.5d seems to demonstrate that there is no difference at the output of the SSPA for ideal and 
non-ideal data using the Butterworth Filter, compared to the plots with no filter and non-ideal data in 
Figures 3.2d and 3.3d where the spurs are obvious. Figures 3.4e/3.4f and 3.5e/3.5f show the output of 
the 4th order Butterworth with ideal and non-ideal data. 

If the amount of attenuation from baseband filtering is examined by comparing with Figures 
3.2c/3.2d and 3.3c/3.3d with 3.4c/3.4d and 3.5c/3.5d (Butterworth Filter), it can be seen that at 
approximately ± 0.5R B , the Butterworth filter attenuates the sidebands by approximately 23 dB. At 
approximately ± 2.5 R B , the attenuation is around 50dB below the peak of the main lobe, placing the 
absolute level of the sidebands under 50 dB after this point. 
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Figure 3.5 - Power Spectra of Filtered 8PSK Baseband Non-Ideal Data 
using 5th Order Butterworth Filter (BT=1) & SSPA at Saturation Level 



Tables 3.3 and 3.4, at the end of this Chapter, summarize the attenuation spectrum levels relative 
to the first data sideband for Ideal and Non-Ideal data. This information about sideband attenuation 
will be important when discussing the utilization ratio (p) and band utilization in a latter section. 


SIMULATION (5) AND (6 > - 3 rd ORDER BESSEL FILTERED BASEBAND DATA 

(IDEAL AND NON-IDEAL DATA) 

Again these two simulations used the same parameters as those discussed for the 5th Order 
Butterworth filter. In the case of the Bessel filter the effect of in-band spurious emissions are more 
present and more evident than for the Butterworth filter as shown in Figure 3.7a compared to Figure 
3.6a. The output of the SSPA (Figures 3.6c/3.6d and 3.7c/3.7d for ideal and non-ideal data 
respectively) show some small spurs emissions caused by the non-ideal data. 

With respect to the sideband attenuation, the values of attenuation for the Bessel (3rd Order) are 
comparable to the 5th Order Butterworth filter (see Tables 3.3 and 3.4). 
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SIMULATION. (7) TO (12) - SRRC FILTER WITH IDEAL AND NON-IDEAL DATA 

The SRRC block in SPW was used as the transmitter and receiver to simulate the SRRC filter as a 

spectrum filter. To simulate a SRRC transmitter filter non-sampled data were used by cascading the 

SRRC with a 1/sinc filter ( function already implemented in SPW). For the receiving SRRC no 1/sinc 

filter was cascaded. Three different rolloff factors were used during the simulation to see the 

difference in bandwidth. The rolloff factors that were used were a = 0.25, 0.5 and 1 . In the SRRC 

filter block, the tap length or FIR length was set to 256 which takes into account 8 intersymbol 

interferences from neighboring pulses since the sampling frequency, f s , was set to 16 samples/sec: 

256 taps = 8 samples interfering before and after. 

2*(fJ 

It is to be noted that the sampling frequency was changed or lowered from 250 samples/sec to 1 6 
samples/sec. Nothing is gained by oversampling (only higher frequencies are seen in the graphs). By 
having a high f s , simulations take more time to terminate. 

Therefore, the following parameters were used for the simulations performed with the SRRC filters: 

- f s = sampling frequency = 16 samples/sec (gives the same results as using 250 samples/sec 
and it is faster) 

- R s = symbol rate = (1/3) symbols/sec; 

- FFT Points = 16384; 

- # samples = 20K; and 

- Tap length (FIR length) = 256. 

Also for data asymmetry a 2% value was used with a probability of zero equal to 0.45. 

Thus in this section, spectra for SRRC with a 1/sinc filter in cascade at the transmitter are shown. 
The SRRC filters have advantages and disadvantages with respect to Butterworth and Bessel filters 
which can only be demonstrated by comparing the several types. Thus, to make the comparison, 
power spectra using several different SRRC filters must be examined. 
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Square R oot Raised Cosine Filter (a = 0.251. 

Comparing Figures 3.8a and 3.9a of a SRRC filter with a roll-off factor of a = 0.25 with Figures 
3.4a and 3.5a for Butterworth and Figures 3.6a and 3.7a for Bessel Filters, it is immediately noted that 
the SRRC filters offers a significant advantage over both of the Butterworth and Bessel filters (BT=1 
for both these filters). The spectrum is narrower, better defined and cleaner than the spectra for other 
filter types. On the other hand, the attenuation is less than that for the Butterworth or Bessel filters at 
high R b (~ ± 20 R b ) as indicated in Tables 3.3 and 3.4. In-band spurious emissions are absent, even 
for non-ideal data. This is the narrowest bandwidth SRRC filter that is considered in this study. 

Outputs from the power amplifier (see Figures 3.8b/3.8c and 3.9b/3.9c for ideal and non-ideal data 
respectively) shows that the nonlinearity (and non constant-envelope signal going in) in the SSPA 
causes an increase in amplitude in the 1 to 4 R B region. Finally Figures 3.8d/3.8e and 3.9d/3.9e shows 
the output of the 4th Order Butterworth Filter. Note how the attenuation remains fairly constant around 
-7R B with a value of approximately -68 dB. 
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Figure 3.8 - Power Spectra of Filtered 8PSK Baseband Ideal Data 
using SRRC (a = 0.25) Filter & SSPA at Saturation Level 





Figure 3 .9 - Power Spectra of Filtered 8PSK Baseband Non-Ideal Data 
using SRRC (a = 0.25) Filter & SSPA at Saturation Level 





Square Root Raised Cosine Filter fra = 0.51 

Similar power spectrum plots are found for the case of a = 0.5 (Figures 3.10a and 3. 1 la) compared 
to Figures 3.8a and 3.9a. The bandwidth is slightly wider for this case (a = 0.5). Nonetheless as with 
SRRC ( a = 0.25 ), the transmitted power spectrum (Figures 3.10e or 3.1 le for ideal and non-ideal 
data respectively) is smooth and the absence of in-band spurious emissions is obvious, even for non- 
ideal data But again, the attenuation produced by the Butterworth or Bessel Filter is much superior at 
higher frequencies than the SRRC (refer to Tables 3.3 and 3.4). Therefore since there appears to be no 
significant difference in the spectra between the SRRC with a = 0.25 and a = 0.5, the selection must 
depend upon other factors such as ISI, implementation complexity, etc. 

Square Root Raised Cosine Filter fa = D 

A SRRC filter with an a = 1 has a larger bandwidth than filters with smaller values of a, therefore 
the baseband data spectrum using ideal and non-ideal data (Figures 3.12a and 3.13a) are wider (than 
for a = 0.25 and a = 0.5). Also, it is to be noted in Figures 3. 12b/3. 12c and 3. 13b/3. 13c, the higher 
order sidebands ( 2nd, 3rd, .. lobes) are more visible in contrast to the other SRRC filters. 

Nonetheless, again at the power amplifier’s output (Figures 3.12b/3.12c or 3.13b/3.13c for ideal 
and non-ideal data respectively), the spectrum is characterized by a smooth roll-off with no-in band 
spurious emissions evident. Due to this filter’s wider bandwidth, the attenuation is slightly less than 
that for filters with smaller values of a especially with the non-ideal data (refer to Tables 3.3 and 3.4). 
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Figure 3.1 1 - Power Spectra of Filtered 8PSK Baseband Non-Ideal Data 
using SRRC (a = 0.5) Filter & SSPA at Saturation Level 




Figure 3.12 - Power Spectra of Filtered 8PSK Baseband Ideal Data 
using SRRC (a = 1) Filter & SSPA at Saturation Level 





Figure 3.13 - Power Spectra of Filtered 8PSK Baseband Non-Ideal Data 
using SRRC (a = 1 ) Filter & SSPA at Saturation Level 





Table 3.3 - Spectrum Levels Relative to First Data Sideband 
(SSPA at Saturation Level and Ideal Data) 


FILTER TYPE 

±or b 

dB 

±5Rb 

dB 

± io R b 

dB 

±20 R b 
dB 

None, Unfiltered Data (Reference) 



-44 

-44 

Butterworth, 5th Order (BT=1) 

0 

-63 


-81 

Bessel, 3rd Order (BT=1) 

0 

-65 

-72 

-80 

Square Root Raised Cosine (a = 0.25), Sampled Data 

0 

-65 

-67 

-71 

Square Root Raised Cosine (a = 0.5), Sampled Data 



-67 

-70 

Square Root Raised Cosine (a = 1), Sampled Data 

0 

-63 

-66 

-70 


Table 3.4 - Spectrum Levels Relative to First Data Sideband 
(SSPA at Saturation Level and Non-Ideal Data) 


FILTER TYPE 
None, Unfiltered Data (Reference) 

Butterworth, 5th Order (BT=1) 

Bessel, 3rd Order (BT=1) 

Square Root Raised Cosine (a - 0.25), Sampled Data 
Square Root Raised Cosine (a = 0.5), Sampled Data 
Square Root Raised Cosine (a = 1 ), Sampled Data 


± 0 R 
dB 
0 


0 


0 


0 


0 


0 


±5 R 
dB 
-39 


-63 


-66 


-67 


-66 


-63 


± 10R 
dB 


-42 


-73 


-74 


-68 


-68 


-67 


± 20 R 
dB 
-43 


-82 


-83 


-72 


-71 


-70 
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POWER SPECTRA FOR NON-FILTERED AND FILTERED g PSfC 
WLTH SSPA BACKED OFF FROM SATURATION 

All of the previous simulations were repeated for the SSPA backed off from saturation by 10 dB, so 
that it is operating in the linear region. Results of these simulations are given in Figures 3.14 to 3.25. 
Since the system is operating in the linear region, the power spectra plots remain the same through the 
SSPA, and the 2nd Harmonic Butterworth filter. 
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Figure 3.17- Power Spectra of Filtered 8PSK Baseband Non-Ideal Data 
using 5th Order Buttenvorth Filter (BT=1) & SSPA at 10 dB Backoff 
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Figure 3.23 - Power Spectra of Filtered 8PSK Baseband Non-Ideal Data 
using SRRC (a= 0.5) Filter & SSPA at 10 dB backoff 
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3.3 End-to-End System Performance: Symbol Error Rate (SER). 

To measure the end-to-end system performance, the symbol error rate was measured for different 
values of average symbol energy to noise ratios (Es/No). Such a measurement was performed using 
SPW and the system described previously. The three types of baseband filters described earlier were 
used. To choose the best filter from such a study one selects the filter that has a BT product such that 
the filter ISI loss is < 0.4 dB [4]. To obtain a measure of the filter IS! loss, a reference needs to be 
established. An ideal linear system is used where ideal data and hardware are utilized to provide a 
reference for comparisons with the system being studied (a plot of Probability of error versus Es/No is 
derived and called Theoretical plot in this case for 8PSK). The proper reference system to measure 
filter loss would be a system where measurements are made by using non-ideal components and data. 

From the experience of the NMSU Center for Space Telemetering and Telecommunications 
Systems, the effect of non-ideal data and system are negligible overall and usually the error produced 
by such non-idealities are approximately 0.5- 1.0 dB. To see how small this effect is, a simulation was 
performed using no baseband filter (BT= oo), the non-linear SSPA and non-ideal data. The SER curve 
is plotted as a function of the SNR. This plot was then compared with the ideal and theoretical curves 
described above and shown in Figure 3.26. The theoretical curve was plotted using the Symbol Error 
Performance expression derived from [8], i.e., for equally likely coherent M-ary PSK signaling: 



where P E (M) is the probability of symbol error, E s = E b (log 2 M) is the energy per symbol, and M=2 X 
is the size of the symbol set i.e X~3 then M=8 since 8PSK signaling is used. 
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From Figure 3.26 it can be seen that the difference between the theoretical curve and the non-ideal 
one is very small. Around SER=10 , the SNR difference between the two curves is approximatelv 0.2 
dB (15.873 dB - 15.68 dB). Note that the last portion of the curve was extrapolated. For this reason 
the simulations for the SER curves were performed using ideal-data (which overlaps with the 
theoretical curve) keeping in mind that there will be a small SNR deviation due to the non-idealities of 
the system. In the CCSDS reports ([3], [4]), simulations were made using non-ideal data and non-ideal 
system components. 

These non-idealities reflect those found in specific equipment and may vary from system-to- 
system. Perhaps measurements of the non-idealities could be made at a later time and then 
incorporated into the overall results of the project. The system performance was then done with ideal 
data (SSPA at saturation level) for the Butterworth and Bessel filter set to BT=1,2 and 3 and also the 
SRRC filters set to a = 0.25, 0.5 and 1 (worse case for SRRC filter). 
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The symbol erTor rate was then plotted as a function of the average symbol SNR. Table 3.5, 
summarizes the simulations performed for the 8PSK SER using the system mentioned above: 

Table 3.5 - SER for 8PSK Simulations 


SIMULATION 

NUMBER 

MODULATION 

TYPE 

FILTER 

LOCATION 

FILTER 

TYPE 



1 

8PSK 

none 

none 

none 

Ideal NRZ-L 

2 

8PSK 

Baseband 

Butterworth 

(BT=1) 

5th Order 

Ideal NRZ-L 

3 

8PSK 

Baseband 

Butterworth 

(BT=2) 

5th Order 

Ideal NRZ-L 

4 

8PSK 

Baseband 

Butterworth 

(BT=3) 

5th Order 

Ideal NRZ-L 

5 

8PSK 

Baseband 

Bessel 

(BT=1) 


Ideal NRZ-L 

6 

8PSK 

Baseband 

Bessel 

(BT=2) 


Ideal NRZ-L 

7 

8PSK 

Baseband 

Bessel 

(BT=3) 

3rd Order 

Ideal NRZ-L 

8 

8PSK 

Baseband 

SRRC 

a = 0.25 

.... 

Ideal NRZ-L 

9 

8PSK 

Baseband 

SRRC 

P 

II 

O 

Lt 

Ideal NRZ-L 

10 

8PSK 

Baseband 

SRRC 

P 

II 

© 

Ideal NRZ-L 


Figures 3.27, 3.28 and 3.29, respectively, show the SER for 8PSK: 5th Order Butterworth BT=1, 2 
and 3 (SSPA at Saturation & Ideal-Data), SER for 8PSK: 3rd Order Bessel BT=1, 2 and 3 (SSPA at 
Saturation & Ideal-Data) and SER for 8PSK: SRRC a=I (no SSPA). a=0.25, 0.5 and 1 (SSPA at 
Saturation & Ideal-Data). 
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Figure 3.29 - SER for 8PSK: SRRC a = 1 (no SSPA), a = 0.25 , 0.5 and 1 
(SSPA at Saturation & Ideal Data) 

From these simulations it is then possible to determine the ISI loss due to the system. This is found 
by taking the difference between the Es/No values measured with and without the filter. As indicated 
in [4], the losses are measured over values of 10' 3 < SER < 10* 2 which is the normal operating region 
for CCSDS encoded data. Losses are then tabulated for BT =1,2,3 (see Table 3.6). Also as indicated 
in that same report the optimal filter is found by selecting the lowest BT providing acceptable ISI loss 
(in this case a threshold of ISI loss < 0.4 dB). The following table gives the results of the baseband 
filter ISI losses at 10 3 SER using the Figures above: 

Table 3.6 - Baseband Filters ISI Loss Measurements at 10° SER 


FILTER TYPE 

LOSS (BT=1) 
(dB) 

LOSS (BT=2) 
(dB) 

LOSS (BT=3) 
(dB) 

Butterworth 5th Order 

17.86-15.68 = 

2.18 

16.294-15.68= 

0.614 

16.048-15.68 = 
0.368 

Bessel 3rd Order 

17.515- 15.68 = 

1.835 

16.9734- 15.68 = 

1.2934 

16.088- 15.68 = 

0.408 

SRRC (a= 1.0) 

16.498 -15.68 = 

0.818 

— 
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Note that only the loss of the SRRC with a=l .0 is shown in the table since the other SRRC filters have 
higher values of losses (as shown in Figure 3.29) and their bandwidths are much smaller than BT=1 
therefore they can not be compared with the Butterworth or Bessel filters. 

Again as stipulated in [4], the optimum filter is the one that has the smallest BT value providing 
losses < 0.4 dB. Therefore with the results shown in the previous table, the following can be stipulated 
as to which filter is the optimum: 

Table 3.7 - Optimum Baseband Filter 


FILTER TYPES 

OPTIMUM BT VALUE 

Butterworth (5th Order) 

BT = 3 

Bessel (3rd Order) 

n BT = 3 

SRRC (a=1.0) 

none 


It must be emphasized that these simulations were performed with ideal data therefore since these 
values of BT barely meet the ISI < 0.4 requirement, with a non-ideal system this threshold would not 
be met (add approximately 0.5 to 1.0 dB to the values measured). Thus a simulation with a higher 
order of BT should be performed (example BT = 4) which would certainly meet the above requirement 
since the bandwidth would be increased. 

In Phase II [3], a Band Utilization Ratio (p) was defined. The interpretation of this ratio can be 
explained as follows. First with the investigation on the baseband filters performed with the 8PSK 
modulation, it w'as noticed that such filters significantly narrow the RF spectrum bandwidth. 

As shown in Tables 3.3 and 3.4, reductions of approximately 24 to 28 dB and 26 to 31 dB were 
achieved at ± 5R B and ± 10 R B respectively. As stated in [3] “these reductions were obtained with 
ideal data (non-ideal data gave similar results) and with no effort devoted to filter optimization. With 
some effort, greater attenuations may be obtainable.” 
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From this, two questions were raised in Phase II about the significant change in the RF spectrum: 

(1) - Can the number of spacecraft using a specific frequency band be increased if 

baseband filtering is used ?; and 

(2) - If the number of spacecraft using a specific frequency band can be increased [ how 
many more can be accommodated than would be the case if no filtering is used? 

To answer to these questions, a frequency band Utilization Ratio (p) can be defined as : 

p = Nymber Q f Spacec raft w i t h Fi ltering Accommod ated in Frequency Band 

Number of Spacecraft without Filtering Accommodated in Frequency Band 

An important assumption in deriving this ratio was made: the “spectra from spacecraft in adjacent 
channels will be permitted to overlap one another provided that, at the frequency where the overlap 
occurs, the signals are at least 50 dB below that of the main telemetry lobe (1st data sideband).” [3]. 

Using SPW, measurements were made using unfiltered and non-ideal data (the 2nd harmonic filter 
was removed) to determine the frequency at which the data spectrum was 50 dB below the main lobe. 
The Power Spectrum plots shown earlier were used to determine this 50 dB level. For example, for 
ideal data without a spectrum shaping filter the -50dB point is situated at approximately 35 R B (refer 
to Figure 3.2d) and using the Butterworth Filter (BT=1) with Ideal Data, the -50 dB point is 
approximately at 2.5 R B (refer to Figure 3.4c), then the Utilization ratio, p, is equal to 


P = 


35*b 

25 R d 


= 14. 


Note that the way the ratio is calculated in this example is different from the formula 
but the result is the same. The number of spacecraft with filtering accommodated 
band will be larger than if it was not filtered. 


given previously 
in the frequency 
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The following table shows the results of Utilization ratio, p, for the other filters. The same 


calculation as performed above were done to complete this table. 

Table 3.8 - Utilization Ratio Improvement for Various Spectrum Shaping 


FILTER TYPE 

Ideal Data 
-50 dB pt. 

Ideal Data 
Util. Ratio 
(P) 

Non-Ideal 

Data 

-50 dB pt. 

Non-Ideal 

Data 

Util. 

Ratio 

(p) 

None, Unfiltered Data (Reference) 

35 R b 

1 

40 R b 

1 

Butterworth, 5th Order (BT=1) 

2.5 R b 

14 

2.3 R B 

17.39 

Bessel, 3rd Order (BT=1) 

2.95 R B 

11.86 

. 

2.5 R b 

16 

Square Root Raised Cosine 
(a = 0.25) Sampled Data 

1.7 R B 


17 R B 

23.5 

Square Root Raised Cosine 
(a = 0.5) Sampled Data 

1.8 R b 

19.4 

1.7 R b 

23.5 

Square Root Raised Cosine 
(a = 1) Sampled Data 

2.5 R B 

14 

2.2 R B 

18.18 


From this last table it is obvious that baseband filtering offers a significant improvement in the 
number of spacecraft operating in a frequency band. This improvement is also superior for the case of 
non-ideal data. In fact by baseband filtering, the bandwidth utilization can increase by a factor of 12 
to 24. 

As mentioned in [3] : “While these ratios may represent upper bounds, it is clear that a significant 
increase in spectrum utilization is potentially possible using baseband filtering, even if the ratios are 
only 3 or 4 to 1 . Judging by the auctions now underway in several countries, such additional frequency 
spectrum will be worth a small fortune and should easily justify expenditures necessary' to develop any 
new filters and/or earth station equipment.” 
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Finally to complement this section on SER, it can be mentioned that different techniques can be 
used to mitigate the effects due to the nonlinear bandlimited signal distortion present in these high data 
rate satellites and therefore improve the End-to-End System Performance (SER). [11] presents a 
development of the low-pass Volterra discrete-time model for a nonlinear communication channel. 
More specifically, the dissertation in [11] is a study of adaptive Volterra equalizers on a downlink- 
limited nonlinear bandlimited satellite channel. It is found that a receiver consisting of a fractionally- 
spaced equalizer (FSE) followed by a Volterra equalizer (FSE- Volterra) gives improvement beyond 
that gained by the Volterra equalizer. 44 Significant probability of error performance improvement is 
found for multilevel modulation schemes. Also, it is found that probability of error improvement is 
more significant for modulation schemes, constant amplitude and multilevel, which require higher 
signal to noise ratios (i.e., higher modulation orders) for reliable operation”. Furthermore, as 
discussed in [12], in systems that have a strict power limitation and high data rate demands (example: 
small satellite transmitting through NASA’s TDRSS network), the dominant noise is present in the 
uplink rather than the downlink channel. Therefore [12] shows that by modifying the standard 
receiver for several digital modulation formats and through proper specification of receiver decision 
regions, pulse shape and receiver filters, significant gains in performance can be produced. On the 
other hand with respect to the bandwidth utilization, [13] introduces a way to develop technology and 
equipment that would allow users to maximize their data throughput in an RF spectrum. All these 
experiments where conducted in the Center for Space Telemetering and Telecommunications Systems 
at New- Mexico State University . 
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3.4 Non-Constant Envelope 8PSK Simulations 

As discussed in the filter section, a finite-duration signal must have infinite bandwidth. Therefore 
when the signal is filtered (to minimize the bandwidth) and the cutoff frequency of the filter is set 
smaller and smaller, it was found with simulations that the ISI increases and the eye diagram 
eventually closes at a point where it exceeds the channel noise. Also, the signal envelope is not 
constant anymore as it will be shown in Figure 3.33. 

The two plots of Figure 3.30 (Figure 3.30a and 3.30b) show the real and imaginary parts 
respectively of the output of the 8PSK modulator. From Figure 3.31a, it is clear that the envelope of 
the signal is constant (flat line with an amplitude of -J2). Figure 3.32 shows, the effect of filtering the 
signal using a Butterworth Filter 5th Order, BT=1. Note how the signal envelope acquires variations 
in Figure 3.33a and how the transitions during each bit are smooth (1 bit corresponds to 16 samples 
since f s = 16 samples/bit). 

Finally if the output signal from the spectrum filter (for example in this case, the 5th Order 
Butterworth Filter), goes through a non-linear device ( a SSPA or TWT for example) , the signal will 
be modified in such a way that the performance of the receiver filter will be affected due to the AM- 
AM and AM-PM conversion produced by the amplifier (this was discussed in the SSPA section). The 
SSPA is a linear amplifier when operated at the appropriate drive level. If the drive level is increased, 
the efficiency (RF output/dc input) is improved, but the amplifier becomes nonlinear. In this case 
constant envelope signals, such as PSK or PM, need to be used so that the intermodulation distortion 
will not cause a problem. This is the mode of operation that is used on the SSPA. Also as mentioned 
in [8], many analysis techniques on non-linear distortion have been published in the literature but none 
of them has been really satisfactory. Figures 3.34 and 3.35 show the output of the SSPA after the 
output of 5th Order Butterworth (BT=1) was used as input. 


73 



Thus when PSK (non-filtered and constant envelope) goes through a high powered amplifier (in the 
non-linear or linear region), the same waveform at the input will be reproduced at the output due to 
the constant envelope of the data. This can be seen in Figures 3.36 and 3.37, where no spectrum 
shaping filter was used thus the signal going through the SSPA has a constant envelope. 
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Figure 3.30a - Output of 8PSK Modulator: Real Part 
Output of 8PSK Modulator: Imaginary Part 
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Figure 3.30b - Output of 8PSK Modulator: Imaginary Part 
Figure 3.30 - Output of 8PSK Modulator: Real & Imaginary Part 



Output of 8PSK Modulator: Real Part 
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Figure 3.31b - Output of 8PSK Modulator: Phase (degrees) 
Figure 3.31 - 8PSK Modulator Output: Magnitude & Phase 
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Figure 3.32b - Output of 5th Order Butterworth Filter (BT=1): Imaginary Part 
Figure 3.32 - 5th Order Butterworth Filter (BT=1) Output: Real & Imaginary' Part 
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Output of 5th Order Butterworth Filter (BT=1): Magnitude 



Figure 3.33a - Output of 5th Order Butterworth Filter (BT=1): Magnitude 
Output of 5th Order Butterworth Filter (BT=1): Phase (degrees) 



Figure 3.33b - Output of 5th Order Butterworth Filter (BT=1): Phase (degrees) 
Figure 3.33 - 5th Order Butterworth Filter (BT=1) Output: Magnitude & Phase 


78 











Figure 3.35a - Output of SSPA (Saturation Level): Magnitude 
Output of SSPA (saturation level): Phase (degrees) 



Figure 3.35a - Output of SSPA (Saturation Level): Phase (degrees) 

Figure 3.35 - Output of SSPA (Input = Output of 5th Order Butterworth Filter (BT=1)): 

Magnitude & Phase 
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Amplitude 
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Output of SSPA (saturation level) No Filter: Real Part 



Figure 3.36a - Output of SSPA (Saturation Level) No Filter: Real Part 


Output of SSPA (saturation level) No Filter: Imaginary Part 



Figure 3.36b - Output of SSPA (Saturation Level) No Filter: Imaginary Part 

Figure 3.36 - Output of SSPA (Input = Output of 8PSK Modulator (Constant Envelope)): 

Real & Imaginary Part 


81 





Figure 3.37a - Output of SSPA (Saturation Level) No Filter: Magnitude 



Figure 3.37b - Output of SSPA (Saturation Level) No Filter: Phase (degrees) 
Figure 3.37 - Output of SSPA (Input = Output of 8PSK Modulator (Constant Envelope)): 

Magnitude & Phase 
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When amplitude modulation is introduced in the input data ( which can be introduced into the PSK 
signal by pulse shaping - filtering), then the output from the SSPA(at saturation level) will be different 
from the input, i.e., if ones tries to amplify an AM signal with this type of SSPA (in the non-linear 
zone) or other nonlinear amplifiers, the AM on the output will be distorted. However, RF signals with 
a constant envelope, such as FM signals, may be amplified without distortion since a nonlinear 
amplifier preserves the zero-crossings of the input signal. Therefore pulse shaping gives a smaller 
bandwidth but a non-constant envelope which reduces the performance of the communication system. 
On the other hand, no pulse shaping gives a larger BW but the signal has a constant envelope which is 
a plus on the performance of the receiver filter. It would then be interesting to investigate the trade-off 
between constant envelope and BW. 

A proposed question was: 

How much variation in the constant envelope occurs 
as the bandwidth is made smaller due to pulse shaping? 

To answer to this question, it is necessary to measure the bandwidth and envelope variation for 
different types of pulse shaping. It would then be possible to see what type of relationship these two 
parameters have. To analyze the effect of non-constant envelope going through the SSPA, different 
simulations where conducted on SPW with the 3 types of filters used for the power spectrum plots and 
SER curves performed in the first and second parts of this report: 

(1) - 5 th Order Butterworth Filter; 

(2) - 3 rd Order Bessel Filter; and 

(3) - SRRC Filters. 
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Figure 3.38 shows a block diagram of the circuit used to perform the simulations (Appendix C 


contains more detailed diagrams and explanations). 



Figure 3.38 - Non-Constant Envelope Simulations Block Diagram 
Note that all simulations were performed at baseband and the SSPA was used at saturation level. 
Figure 2.10 showed the magnitude characteristic curve of the ESA 10 Watts SSPA that was used in 
these simulations. Figure 3.39 shows the same curve as Figure 2.10 but zoomed in the non-linear and 
saturation level regions (between -6 to 1 dB for the Input Level Power). 



Figure 3.39- SSPA, ESA 10 Watts, Zoomed Magnitude Characteristic Curve 
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Note that that the saturation level region and non-linear region are very close, i.e., approximately -2 dB 
which corresponds to a voltage of 0.794 ( 10 ( ' 2dB/20) ). A voltage of 1 is at OdB input level Power. 
Therefore all simulations were performed using a reference of 0 dB or 1 volt. 

In Figure 3.40, the Phase Response of the SSPA that was used is shown. 



Appendix B contains the table of the values of this amplifier and the block diagram used to 
simulate or interpolate the values from the table. From the phase response curve (Figure 3.40), it is to 
be noted that a phase rotation of approximately -5.4° occurs at 0 dB of input level. This will cause a 
rotation of the 8PSK. symbols as shown in Table 3.9. 
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Table 3.9 - Rotation of 8PSK Symbols due to PM Conversion on SSPA 


p 

M CONVI 

iRSlON AT SATURATION LEVEL (~0dB Input Level) 

SYMBOLS 

1 

2 

3 

4 

5 

6 

7 

8 

Initial Position (°) 

22.5° 

67.5° 

112.5° 

157.7° 

202.5° 

247.5° 

292.5° 

337.5° 

PM Rotation (°) 

n 

H 

-5.4° 

b 

H 

B 

B 

B 

Final Position (°) 
= Initial Position 
- PM Rotation 

I7.1.° 

62.1° 

107.1° 

— 

152.3° 

197.1° 

or 

-162.9° 

242.1° 

or 

-117.9° 

287.1° 

or 

-72.9° 

332.1° 

or 

-27.9° 


From Table 3.9, a rotation of -5.4° will be introduced to any signal with 0 dB gain (on a 1 volt 
reference). Thus if a non-constant envelope signal goes through this SSPA, a maximum rotation of 
5.4° will be induced but also smaller values of rotation will be introduced ( < | -5.4° |) depending on 
the input level of the non-constant envelope signal. For example, Figure 3.41 shows the scatter plot 
of Figures 3.30, 3.32, 3.34, i.e., the 8PSK Modulator, 5th Order Butterworth Filter (BT=1) and SSPA 
(Saturation Level) Outputs. A scatter plot, also known as the signal “constellation” or “I-Q plot”, 
shows the location in the complex plane of the in-phase and quadrature-phase components of a 
complex signal at discrete time samples. 
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Scatter Plot: 8PSK Modulator, 5th Order Butterworth Filter (BT=1) and SSPA Output (saturation level) 



real 

Figure 3.41 - Scatter Plot: 8PSK Modulator, 5th Order Butterworth Filter (BT=1) 
and SSPA (Saturation Level) Outputs 

Note that in this last figure the symbol “x” represents the symbols output from the 8PSK Modulator. 
The position of these symbols is exactly at [(7i/8)*I + (ti/ 4)] and have an amplitude of 1 (everything 
was normalized to 1 volt or OdB). The second set of points (cluster of points outside the unit circle) 
represents the output of the 5th Order Butterworth Filter with BT=1. Note how the symbols are not 
shown as a single element anymore; this is due to the non-constant envelope introduced by filtering 
the input signal (output of the 8PSK Modulator). This cluster will be larger or smaller depending on if 
the signal has a more constant envelope or non-constant envelope, respectively. The signal was not 
rotated or phase shifted going through the filter therefore the mean of the cluster is approximately 
[(n/8)*i + (ti/4)]. Finally the last set of points (points aligned on the unit circle) represents the output 
of the SSPA. 
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Note that this time the symbols have been rotated (* -5.4°) due to the AM/PM conversion 
introduced by the SSPA. Also the symbols are aligned on the unit circle since the SSPA is operating at 
saturation level (0 dB or 1 volt). 

To be able to create this scatter plot from Figures 3.30, 3.32 and 3.34, the following procedure was 
followed. First a simulation on SPW was performed using the block diagram shown in Figure 3.34 (a 
more detailed diagram in shown and explained in Appendix C). Then when the simulation was 
completed, the American Information Code for Information Interchange (ASCII) file containing the 
complex values of the output of the SSPA for example (signal shown in Figure 3.30) was transferred 
(File Transfer Protocol (FTP)) to a Personal Computer (PC 486 66 MHz) where the signal was 
sampled to its optimum value to produce the SSPA Output scatter plot shown in Figure 3.41. For 
example, for the signal in Figure 3.30 , a signal delay equal to 20 samples had to be used to get the 
optimum scatter plot. This delay of 20 samples was calculated as follows: 

First the pulse shaping filter (5th Order Butterworth Filter (BT=1)) and SSPA introduced a delay of 
12 samples (includes the delay due to the transition) and since the sampling of the received bit has to 
occur in the middle of the bit ( f s /2 where f 5 = sampling frequency = #samples/bit), then 

T sampling time = (Delay of 5th Order Butterworth filter (BT=1) and SSPA) + (f s /2) 

where for the 5th Order Butterworth Filter (BT=1), the sampling frequency = f s = 16 
samples/bit (for more information on the choice of the value of f s , refer to Appendix C). 

Thus 

T sampling time = 1 2 samples + (1 6 samples/bit 12 ) = 12 + 8 


T sampling time ^0 Samples 
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All this was processed in a MATLAB n file called SPWSOBRE.M (see Appendix D). This 
procedure reads the ASCII file transferred from SPW and samples it as discussed above to produce a 
scatter plot. Afterwards, the mean and variance of each symbol is calculated. 

The mean for each symbol, in one of the eight decision regions (described below) is calculated as 
follows: 


where 


mean{i) 


= mean[Q(i )] 


for i = 0 to 7 inclusively 


•0(0 are all the angles of the 
delimited by: 

region 0: 
region 1: 
region 2: 
region 3: 
region 4: 
region 5: 
region 6: -ix/2 <0 
region 7: -rc/4 <0 


71/4 

or 

0 <0 

<45°; 

n/2 

or 

45° <0 

<90°; 

6ti/8 

or 

90° <0 

< 135°; 

71 

or 

135° <0 

< 180°; 

-6tt/8 

or 

CD 

V 

o 

O 

oo 

I 

<-135°; 

-71/2 

or 

-135° <0 

< -90°; 

-tt/4 

or 

-90° <0 

< -45°; and 

0 

or 

-45° <0 

<0°; 


symbol contained in one of the eight decision regions 

0 < 0 
7l/4 <0 
n/2 <0 
6tt/8 <0 
-71 <0 

-6tt/8 <0 


+t MATLAB is a registered trademark of The MatWorks Inc., Cochituate Place, 24 Prime Park Way, 
Natick, Mass. 01760. 
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• mean[x] is a function in Matlab that calculates the mean of the vector x. In this case x is a 
vector containing the angles of the symbols in one of the 8 decision regions. The mean in 
Matlab is calculated in the following manner: 


1 

mean[Q(i)] = — 


N 




where i = decision region (0 < i < 7); 


x = element in the vector 0 (i) ; and 


N - number of points in the decision region or vector© (i). 


Also the variance of the symbol vectors are calculated to see the distribution of the symbols around 
the mean in the decision region. In SPWSOBRE.M the following line calculates the variance: 

- [sfcf[9(01] 2 ~ 

where std = standard deviation which is a function defined in Matlab as being the 
sample standard deviation normalized by the square root of (N-l), where N is the sequence or vector 
length. Thus in other words, the variance is calculated as 



where 


i = 0 to 7 (decision regions); and 
N = number of points in the decision region. 


90 






Figure 3.42 graphically shows the mean and variance in a decision region: 



Figure 3.42 - Mean and Variance in Decision Region 
Thus, for the 5th Order Butterworth filter with BT=1, the scatter plot due to the SSPA shown in 
Figure 3.41 has the following values of mean and variance (calculated as shown above) for the 8 
symbols: 

Table 3.10 - Mean and Variance of SSPA Output (Saturation Level) 
with 5 th Order Butterworth (BT=1) 


SYMBOLS 

1 

2 

3 

4 

5 

6 

7 

8 

©mean O 

17.020° 

62.250° 

106.80° 

152.58° 

-162.98° 

or 

197.02° 

-117.86° 

or 

242.14° 

-72.650° 

or 

287.35° 

-27.753° 

or 

332.247° 

rj 

u variance 

(xlO‘ 2 )° 

7.1477 

6.3547 

5.2001 

6.6354 

5.4146 

5.4370 

6.6427 

7.4719 
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If the 2 nd row of Table 3.10 is compared with the last row of Table 3.9, it is noted that the results 
are very similar. The same simulations as described above were conducted for different Bandwidth- 
Time (BT) products for the 5th Order Butterworth Filter and the 3rd Order Bessel Filter. 

Also, the SRRC filter was simulated for roll-off factors varying from 0. 1 to 1 (a = 0.0 could not be 
simulated due to software capabilities; the roll-off factor a has to be > 0.0 in SPW). Table 3.11 
summarizes the different bandwidth and spectrum filters that were used to see the effect of non- 
constant envelope through an SSPA. 


Table 3.1 1 - Non-Constant Envelope Simulations: Filters and BW used 


SPEC TRI M FU. H R 

BANDWIDTH 

5 th Order Butterworth 

BT = 1, 1.05, 1.1, 1.2, 1.4, 1.6, 1.8,2,22, 
2.4, 2.6, 2.8, 3, 3.2, 3.4, 3.6, 3.8, 4 

3rd Order Bessel 

BT = 1, 1.05, 1.1, 1.2, 1.4, 1.6, 1.8, 2, 2.2, 
2.4, 2.6, 2.8, 3, 3.2, 3.4, 3.6, 3.8, 4 

SRRC 

a = 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 
0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 
0.9, 0.95, 1.00 


The following figures, i.e., Figure 3.43 - 5th Order Butterworth Filter: Average Symbol Variance vs 
BT, Figure 3.45 - 3rd Order Bessel Filter: Average Symbol Variance vs BT and Figure 3.47 - SRRC 
Filter: Average Symbol Variance vs Roll-off factor (a), show the average symbol variance plotted 
against the bandwidth. The Average Symbol Variance was calculated as follows: 

_ 8 

^ ® var iance(i) 



8 

which can be explained as being the average of the sum of the variances in the decision regions. 
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These figures show the effect of varying the non-constant envelope of a signal ( by using different 
BT or a) through an SSPA. The Average Symbol Variance is used as a parameter indicating the 
spreading of the symbol with respect to the mean for each of the 8 decision regions. For the 5th Order 
Butterworth Filter in Figure 3.43, it is noted that the variance substantially decreases from BT=1 to 
BT=1.1. This can be explained by the fact that for BT=1, the filter is narrower and contains only the 
main lobe of the power spectral density (refer to section 3.2). 

Thus the signal in the time domain will be more distorted (non-constant envelope) than when a 
larger order of BT (or bandwidth) is utilized (the envelope becomes more constant). 


5th Order Butterworth Filter Symbol Average Variance vs BT 



Figure 3.43 - 5th Order Butterworth Filter: Average Symbol Variance vs BT 
It is interesting to note how the variance ripples down as the BT is increased. Eventually as BT is 
very large (at the limit BT= »), the variance should be equal to 0 since almost all of the PSD (or 
energy) of the signal is captured. 
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Figure 3.44 shows a zoomed version, between BT=1.2 and 3.2, of Figure 3.43. Note the symmetry of 
the ripples as they go to 0. The simulations points “o” were linked together using an interpolation 
function called “spline” in Matlab. 


x -iq' 5 5th Order Butterworth Fitter Symbol Average Variance vs BT 



Figure 3.44 - Zoom of Average Symbol Variance vs BT for 5 th Order Butterworth Filter 
Appendix E contains the scatter plots of the points for Figures 3.43 and 3.44; note that one point in 
Figure 3.43 corresponds to the average variance of the symbols in one scatter plot. Also a listing of 
the program (BUTTPLOT.M) written in Matlab, that was used to create the plot, is included. 

Figure 3.45, shows the results of variance vs BT for the 3rd Order Bessel Filter. Again as in the 
Butterworth case presented earlier, the average variance gradually decreases in a sinusoidal form. The 
damping is more pronounced for the 3rd Order Bessel than the 5th Order Butterworth since the 
Bessel’s variance starts to become a flat line around BT=1.8 but for the Butterworth Filter case the 
steady line becomes present around BT=3. Also note that the average variance is appreciably less in 
magnitude for the Bessel filter (order of magnitude * lO' 4 ) than the Butterworth filter (order of 
magnitudes 10’ 2 . Figure 3.46 shows a zoomed version (between BT = 1.4 and 3.7) of Figure 3.45. 
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Variance in Degrees 


x iQ** 3rd Order Bessel Filter Symbol Average Variance vs BT 



Figure 3.45 - 3rd Order Besse! Filter: Average Symbol Variance vs BT 


3rd Order Bessel Filter Symbol Average Variance vs BT 
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Appendix F contains the scatter plots of the points for Figure 3.45 and 3.46; again, note that one 


point in Figure 3.45 or 3.46 corresponds to the average variance of the symbols in one scatter plot. 


Also a listing of the program (BESSPLOT.M) written in Matlab, that was used to create the plot, is 


included. If the scatter plots in Appendix E and F are compared, it is noted that in fact the average 


variance is significantly less for the 3rd Order Bessel filter than the 5th Order Butterworth Filter. 


Finally Figure 3.47 shows the average variance versus the roll-of factor for the SRRC filter. 


SRRC Filter Average Symbol Variance vs Roll-of Factor (alpha) 
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Figure 3.47 - SRRC: Average Symbol Variance vs Roll-off factor (a) 


In this case note that the aspect of the plot is really different from the two previous filters. In fact, 


for the SRRC, the plot appears to have the aspect of a second order parabola with a maximum at 


approximately a=0.45. Note that simulations could not be performed for 0 < a < 0.1 due to SPW 


software limitations. Nonetheless, it can be seen that the curve gradually decreases to 0 as a reaches 


0. This can be explained by the fact that for a=0, the frequency response of a SRRC is a pure square 


wave (or “brick wall”). Thus in the time domain, it would then correspond to a pure sine (sin(x)/x) 
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function which would have the zero crossings at the sampling frequency therefore eliminating any ISI. 
For the case where a reaches 1, the average variance also decreases toward 0 (but it does not reach 
it). This can be explained since the SRRC filter has a bigger BW for higher a therefore allowing 
more of the power spectrum to be included in the transmission; this is the same result as increasing the 
BT for the Butterworth or Bessel Filter. Finally note how the magnitude of the average variance is 
larger for the SRRC filter than the Butterworth and Bessel. 

Appendix G contains the scatter plots of the points for Figure 3.47. Also a listing of the program 
(SRRCPLOT.M) written in Matlab, that was used to create the plot, is included. If the scatter plots in 
Appendix E, F and G are compared, it is noted that in fact the average variance is significantly bigger 
for the SRRC filter than the 3rd Order Bessel filter or the 5th Order Butterworth Filter. 
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Chapter 4 


CONCLUSIONS AND RECOMMENDATIONS 

From the simulations on 8PSK Baseband Filtering performed in the Center for Space Telemetering 
and Telecommunications Systems at NMSU, the following conclusions can be made on the PSD, SER 
and Non-Constant Envelope: 

With respect to PSD and SER for the 8PSK signal which is an extension of the work performed during 
Phase I, II and III (refer to [1], [2], [3] and [4]) of this study : 

• In-band spurious emissions are more present and more evident for the Bessel Filter than the 
Butterworth Filter (BT=1). With respect to the sideband attenuation, it was found that the values 
of attenuation for the 3rd Order Bessel are comparable to the 5th Order Butterworth (refer to 
Tables 3.3 and 3.4 in the report). 

• For SRRC filters with a = 0.25 and a = 0.5, the bandwidth is narrower than the Butterworth and 
Bessel Filters but the attenuation is less at high frequencies. Nonetheless, the absence of spurious 
emissions is a net advantage. For SRRC a = 1, the bandwidth is wider than the two previous 
SRRC filters and the absence of in-band spurious emissions was again noticed. Less sideband 
attenuation was recorded for this roll-off factor compared with a = 0.25 and a = 0.5. 

• For SER, it was found that the Butterworth and Bessel Filters just barely meet the threshold of ISI 
loss < 0.4 dB at SER = 10 ° . Also the SRRC filters do not meet this specification. Since this 
threshold is barely met by the first two filters and this simulation was performed with only ideal 
data and system it seems reasonable that a simulation with BT = 4 should be performed which 
would certainly meet the threshold. 
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From these results, it can be stipulated that the 5th Order Butterworth has a small advantage 
compared to the other filters. Also more simulations should be performed on the power spectrum 
content depending on the value of BT that meets the threshold mentioned above. Therefore, for 
example for BT = 4 , the power spectrum plots should be produced to be able to calculate the new 
Band Utilization ratio, p. 

Overall, it was shown by using baseband filtering the bandwidth utilization can be improved by a 
factor of approximately 12 to 24 with BT=1 and 8PSK (see Table 3.8 in this report) which can 
significantly increase the spectrum utilization. In Phase II (refer to [3]), it was demonstrated that for 
BPSK, this bandwidth utilization can be potentially increased from 6 to 10 times with BT=1. Thus 
with respect to bandwidth, 8PSK gives better results (using the same bit rate). 

To increase our knowledge with respect to a non-linear satellite channel, specific measurements were 
performed on the Non-Constant Envelope of the signals going through the non-linear SSPA. These 
simulations have not been previously reported. 

• The tradeoffs between constant envelope and bandwidth were observed with simulations for the 
5th Order Butterworth Filter, 3rd Order Bessel Filter and SRRC Filter. The more shaped the pulse 
(to decrease the bandwidth), the more non-constant the envelope of the signal becomes therefore 
increasing errors due to non-linearities in the system. 

=> It was found that the amplitude of the Average Symbol Variance of the signal 
using the 3rd Order Bessel is much less than for the 5th Order Butterworth 
Filter. Also the amplitude for the 3rd Order Bessel Filter damps toward zero 
around BT=1.8 where for the 5th Order Butterworth Filter this damping will 
occur around BT=3. These curves have the aspect of a sinusoidal function 
decreasing exponentially toward zero; and 
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=> The Average Symbol Variance vs a for the SRRC was found to be quite 
different than the two previous spectrum shaping filters. This plot looks more 
like a second order parabola with a maximum at approximately a = 0.45. Also 
the amplitude is significantly larger than the two previous filters. 

Suggestions for further study 

• A non-ideal system can be used for future simulations for SER . This non-ideal system should be 
predetermined with respect to what kind of components have to be used and the amount of error 
they introduce. 

• Increase the BT to 4 and perform more simulations for PSD and SER. 

• For the non-constant envelope simulations more simulations should be performed using different 
orders of filters for Butterworth and Bessel: 

=> From these simulations it might be possible to find a general expression on how 
the variance of the symbols is altered by changing BT. This expression should 
have the following form: a*exp(-x)*sin(x) for the Butterworth and Bessel 
Filters; and 

It would then be possible to predict, with a general equation, when the BT has 
no real effect on the shaping of a pulse. 
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APPENDIX A 


SPW DETAIL DIAGRAMS 
FOR 

PSD AND SER SIMULATIONS 

This appendix contains the detailed diagrams of the overall system used for the PSD and the SER 
simulations. Explanations on how these simulations were performed will be given. 

A.1 PSD Detailed Diagram 

Figure A.l shows a detailed block diagram of the system used to simulate the PSD with different 
spectrum shaping filters. First, a MPSK Source block is used and since the simulations were 
performed for 8PSK the corresponding parameters were given (see right top comer of Figure A. 1 ) to 
have an 8 level output. The output of this block is saved in a library “rubensigs/modout” for later use. 
For a more detailed aspect of this block refer to Figure 2.3 in the report. The output of this modulator 
is then introduced in the “Baseband Spectrum Shaping Filter” (for this case the Butterworth Filter is 
shown). Figure A.2 shows the contents of this Butterworth Filter Block. When performing 
simulations the filter order was set to 5 for the Butterworth filter, the attenuation at the passband edge 
was set to 3dB. These two parameters were kept the same for all the simulations. The sampling 
frequency was set to 250 samples/bit to be able to have as much spectrum as possible giving a range 
from -125 R b to 125 R B . 


103 



PARAMETERS MPSK MODULATOR PARAMETERS 



Figure A.l - PSD Detailed Diagram 
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Figure A.2 - Buttervvorth Filter Block Diagram 
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For the “Passband Edge Frequency “ parameter, B, the value was set to 1/3 since for this simulation 
the BT of the filter was 1 thus 

BT = B* (Symbol Period) = B*(l/ Symbol Rate) = B/ R s 
= B/ (R b / X) = B*MR b 
where R s = Symbol Rate (symbols/sec); 

R b = Bit Rate = 1 bit/sec (set by default for these simulations); 

X = log 2 8 = 3; and 

B = Passband Edge Frequency or Single-Sided 3dB Bandwidth . 

Thus BT = B*(3/l) = 3B =1 which implies that the Passband Edge Frequency = B = 1/3 as shown 
in Figure A.2 and explained above. If the simulation had to be performed for BT=2, then the Passband 
Edge Frequency, B, would be equal to 2/3 and for BT=3, B=1 . The output of this filter is also saved in 
a file for later use. The next block of Figure A.l is the SSPA which is explained in detail in Appendix 
B. It is to be noted that the output of the SSPA is also saved in a file. Finally in the study of the PSD 
of 8 PSK, a 2nd Harmonic Filter (4th Order Butterworth Filter) follows the power amplifier. This type 
of filter is frequently used after the power amplifier to protect users of other bands. The parameters of 
this filter can be set in the same manner as explained above for the 5th Order Butterworth. Also for 
this filter, the output is sent into a file. 

In the case where non-ideal data had to be simulated in the system the parameter “zeta” or data 
asymmetry in the top left comer box in Figure A.l had to be set in following manner. First the 
parameter 4 Symbol Rate (Hz) or RsynT was set to 1/3 (R B / X) as explained above. Using this result, 
the parameter “N (#samples/symbol = f/Rsym) “ is equal to 






250sampies / sec 
(1/ 3)symbol / sec 


= 750 samples f symbol 
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and the parameter “zeta” is defined as N*0.02 since, as explained in this report, CCSDS limits data 
asymmetry (or zeta in our case) to ± 2% (ratio of time duration of a 1 to time duration of a 0). Thus 
using the formula and result for N shown above: 

Nx0.02 = 750x0.02 = 15 samples/symbol of data asymmetry. 

This number would then be entered next to the zeta parameter. Also another parameter that has to 
be changed when simulating with non-ideal data are the “Probability of Zero” or data imbalance which 
was defined in this report as the probability of a bit 1 versus the probability of a bit 0. CCSDS limits 
data imbalance to 0.45. Thus this value of 0.45 would then be entered for the “Probability of Zero”. 

The “Carrier Frequency Wc (Hz)” parameter was set to 0 for all the simulations since these were 
performed at baseband. Also the “SSPA Backoff (dB)” parameter was set to 0 since the simulations 
were performed at saturation level (except when the 10 dB Backoff simulations were performed; then 
the value of SSPA Backoff was set to 10 dB). 

When these parameters are set, the simulation is started by using the “run” function in SPW. For 
the PSD simulations, 300 000 samples were used. After the simulation is terminated, the FFT function 
in SPW is used (FFT Points was set to 524288 for better accuracy as discussed in section 3.2 of this 
report) to plot the PSD diagrams. The final diagrams of the PSD are shown in this report, section 3.2. 
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A.2 SER Detailed Diagram 


The detailed diagram from SPW that was used to measure the SER and produce the corresponding 
SER Curves is shown in Figure A.3. Basically this circuit takes the output signal of the matched filter 
(since NRZ-L is the data format, a sliding integrator was used as the matched filter. For more 
information about the matched filter refer to section 2.7 in this report) which is the signal that was 
distorted through the channel and compares it with the input signal. This comparison is performed by 
a delay & phase meter which tries to correlate the non-distorted signal with the distorted one. When 
the best correlation is found, the distorted and non-distorted PSK signals are converted to their 
corresponding symbol numbers. Afterwards these symbol numbers are compared in a “Simple Error 
Rate Estimator” to find the Symbol Error Rate (SER). 

All this is performed in real time and is shown on the screen of the Sparc- 10 Sun Station as it can 
be seen in the two boxes “ Display and Control” and “Error Rate Display” in Figure A.3. Finally 
Figure A.4 shows a detailed block diagram for the AWGN that was used in the SPW simulations. 
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Figure A. 3 - SER Detailed Block Diagram 
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SIMULATION PARAMETERS 















Figure A.4 - Detailed Block Diagram for the AWGN 
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APPENDIX B 


SSPA, ESA 10 Watts 

This appendix shows the SSPA model used in SPW. The magnitude and phase plots of the SSPA 
used in the simulations performed with 8PSK modulation with different spectrum shaping filters are 
shown in Figure B. 1 . 


SSPA, ESA 10 Watts Amplifier: Magnitude Response Curve 



SSPA, ESA 10 Watts Amplifier: Phase Response Curve 



Figure B.l - Magnitude and Phase of SSPA, ESA 10 Watts 





The data used to produce these plots is listed in page 1 14 and was obtained from JPL (the same 
amplifier was used by JPL for their simulations at saturation level). The data file consists of three 
columns: the first column is the input power in dB, the second column is the output power in dB, and 
the third column is the output phase in degrees (°). Therefore, the data are formatted in such a way 
that SPW can have access to the different numbers in the file. The block “TWT” in the 
“Communication Library” in SPW was utilized. Figure B.2, shows the aspect of this block. 



Figure B.2 - SPW Amplifier Model: Block Diagram 
This block implements the AM-to-AM and AM-to-PM characteristics for a typical TWT or in this 
case a SSPA. Parameters to this block include the operating point in dB ( 0 dB for saturation level), 
the name of the TWT or SSPA file, and the number of interpolated points ( for the simulation 
performed, 240 interpolated points were used). Figure B.3 - Detailed SPW TWT Model, shows a 
detailed version of Figure B.2. SSPAs or TWTs are generally modeled using AM/AM and AM/PM 
conversion characteristics. The characteristics give output power versus input power (AM/AM 
conversion) and output phase shift versus input power (AM/PM conversion). In this block, these 
characteristics are specified by a file of measured input/output data as mentioned earlier. This data file 
represents, for these simulations, the specifications of the SSPA, 10 Watt, S-band power amplifier 
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TUT ..TABLE BLOCK PARAMETERS 



Figure B.3 - Detailed SPW TWT Model 
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used by the European Space Agency (ESA). As stipulated in [14], “ the measured data in the data file, 
consists of relative power levels, referenced to saturation point. That is, the input power data are the 
input power relative to that which causes saturation at the output, and the output power data are the 
output power relative to the power at saturation. This allows the characteristics to be modeled around 
the point where an input envelope value of unity gives an output envelope value of unity. Thus, in the 
measured data , 0 dB input power corresponds to 0 dB output power.” 



SSPA3 TABLE 


$ 

SIGNAL TYPE = double 
VECTOR TYPE = interlaced 
VECTOR LENGTH = 210 
NUMBER OF VECTORS = 3 
NUMBER OF SIGNAL POINTS = 210 
SAMPLING FREQUENCY = 1 
STARTING TIME = 0 
$ 


-3.0000000e+01 

-2.5340000e+01 

-1.5886600e-02 

-2.9500000e+01 

-2.4840000e+01 

-1.6827900e-02 

-2.9000000e+01 

-2.4340000e+01 

-1.7825000e-02 

-2.8500000e+01 

-2.3840000e+01 

-1.8881200e-02 

-2.8000000e+01 

-2.3340000e+01 

-2.0000000e-02 

-2.7500000e+01 

-2.2840000e+01 

-2.1 185100e-02 

-2.7000000e-K)l 

-2.2340000e+01 

-2.2440400e-02 

-2.6500000e+01 

-2.1840000e+01 

-2.3770000e-02 

-2.6000000e+0I 

-2.1340000e+01 

-2.5178500e-02 

-2.5500000e+01 

-2.0840000e+01 

-2.6670400e-02 

-2.5000000e+01 

-2.0340000e+01 

-2.8250700e-02 

-2.4500000e+01 

-1.9840000e+01 

-2.9924700e-02 

-2.4000000e+01 

-1.9340000e+01 

-3.1697900e-02 

-2.3500000e+01 

-1.8840000e+01 

-3.3576 100e-02 

-2.3000000e+01 

-1.8340000e+01 

-3.5565000e-02 

-2.2500000e+01 

-1.7840000e+01 

-3.7673000e-02 

-2.2000000e+01 

-1.7340000e+01 

-3.9905200e-02 

-2.1500000e+01 

-1.6840000e+01 

-4.2269800e-02 

-2.1000000e+01 

-1.6340000e+01 

-4.4774400e-02 

-2.0500000e+01 

-1.5840000e+01 

-4.7427500e-02 

-2.0000000e+01 

-1.5340000e+01 

-5.0237700e-02 

-1.9500000e+01 

-1.4840000e+01 

-5.3214500e-02 

-1.9000000e+01 

-1.4340000e+01 

-5.6367700e-02 
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-1.8500000e+01 
-1.8000000e+01 
-1.7500000e+01 
-1.7000000e+01 
-1.6500000e+01 
-1.6000000e+01 
-1.5500000e+01 
-1.5000000e+01 
-1.4500000e+01 
- 1 .4000000e+0 1 
-1.3500000e+01 
-1.3000000e+01 
-1.2500000e+01 
-1.2000000e+01 
-1.1500000e-K)l 
-1.1 000000e+0 1 
-1.0500000e+01 
-1.0000000e+01 
-9.5000000e+00 
-9.0000000e+00 
-8.5000000e+00 
-8.0000000e+00 
-7.5000000e+00 
-7.0000000e+00 
-6.5000000e+00 
-6.0000000e+00 
-5.5000000e+00 
-5.0000000e+00 
-4.5000000e+00 
-4.0000000e+00 
-3.5000000e+00 
-3.0000000e+00 
-2.5000000e+00 
-2.0000000e+00 


-1.3840000e+01 

-1.3340000e+01 

-1.2840000e+01 

-1.2340000e+01 

-1.1840000e+01 

-1.1340000e+01 

-1.0840000e+01 

-1.0340000e+01 

-9.8400000e+00 

-9.3400000e+00 

-8.8400000e+00 

-8.3400000e+00 

-7.8400000e-f-00 

-7.3400000e+00 

-6.8400000e+00 

-6.3400000e+00 

-5.8400000e+00 

-5.3400000e+00 

-4.8400000e+00 

-4.3400000e+00 

-3.8400000e+00 

-3.3400000e+00 

-2.8400000e+00 

-2.3400000e+00 

-l.8400000e+00 

-1.3400000e+00 

-8.4000000e-01 

-5.9000000e-01 

-3.7999900e-01 

-2.6000000e-01 

-2.0000 100e-01 

-1.6000000e-01 

-1.2000200e-01 

-7.9999800e-02 


-5.9707600e-02 

-6.3245600e-02 

-6.6993 100e-02 

-7.0962700e-02 

-7.5167500e-02 

-7.962 1400e-02 

-8.4339300e-02 

-8.9336700e-02 

-9.4630200e-02 

-1.0023700e-01 

-1.0617700e-01 

-1.1 246800e-0 1 

-1.1913200e-01 

- 1 .26 191 00e-0 1 

-1.3366900e-01 

-1.4158900e-01 

-1.4997900e-01 

-1.5886600e-01 

-1 .6827900e-01 

-1.7825000e-01 

-1.8812000e-01 

-2.0000000e-01 

-4.0000000e-01 

-6.0000000e-01 

-8.0000000e-01 

-1.0000000e+00 

-1.2000000e+00 

-1.4000000e+00 

-1.6000000e+00 

-1.7500000e+00 

-2.2500000e+00 

-2.7500000e+00 

-3.5000000e+00 

-4.0000000e+00 
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-1.5000000e+00 

-5-OOOOOOOe-Ol 

0.0000000e+00 

5.0000000e-01 

1 .0000000e+00 

1.5000000e+00 

2.0000000e+00 

2.5000000e+00 

3.0000000e+00 

3.5000000e+00 

4.0000000e+00 

4.5000000e+00 

5.0000000e+00 


-4.0001400e-02 
0.0000000e+00 
-5.1771900e-07 
'5.1771900e-07 
-5. 1771900e-07 
-5. 1771900e-07 
-5.1771900e-07 
-5.1771900e-07 
-5.1771900e-07 
-5.1771900e-07 
-5.1771900e-07 
-5.1771900e-07 
-5.1771900e-07 


-4.3500000e+00 

-5.0000000e+00 

-5.4000000e+00 

-5.4000000e+00 

-5.4000000e+00 

-5.4000000e+00 

-5.4000000e+00 

-5.4000000e+00 

-5.4000000e+00 

-5.4000000e+00 

-5.4000000e+00 

-5.4000000e+00 

-5.4000000e-H)0 
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APPENDIX C 


NON-CONSTANT ENVELOPE 
SIMULATION 
BLOCK DIAGRAM 

This appendix contains a detailed diagram of the non-constant envelope simulation circuit used on 
SPW to be able to measure the effect of the non-constant envelope through the SSPA. A brief 
explanation of this circuit will be given. 

Figure C.l shows the detailed SPW diagram of the Non-Constant Envelope Simulations also 
described in Figure 3.34 of section 3.4 in this report. Note the similarity of Figure C.l with the 
Figures in Appendix A. One of the only differences is the model of the SSPA which has been slightly 
modified. In fact. Figure C.2 shows a detailed diagram of this new model of SSPA. This new model 
was created to get a better rotation accuracy due to the AM/PM conversion in the SSPA. In this new 
SSPA model, the SSPA table shown in Appendix B has been modified in the sense that two new tables 
were created ,i.e., one table which was named “outdb.ascsig” contained the input dB column (first 
column of the table in Appendix B) and the output dB (second column of table in Appendix B). 
Therefore this first file was the AM/AM conversion created by the SSPA. The second file that was 
created was called “outphase.ascsig” and contains the columns from the table in Appendix B that 
produced the AM/PM conversion ,i.e., column 1 and 3. 

Therefore the two TABULAR NON LINEAR blocks in Figure C.2 refer to these files to output the 
corresponding amplitude or phase conversion. The highest TABULAR NON LINEAR block in the 
figure refers to the “outdb.ascsig” table and the second block refers to the “outphase.ascsig” table. 
The rest of die block diagram is self explanatory. 
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MPSK MODULATOR PARAMETERS 



Figure C.l - Detailed SPW Diagram of Non-Constant Envelope Simulations 
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SPECTRUM SHAPING 










Figure C.2 - Detailed SPW Diagram Of New SSPA 
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The most important point that was considered during these non-constant envelope simulations using 
the diagram shown in Figure C.l was the selection of the sampling frequency parameter, f SJ first 
parameter in the top left comer of Figure C.l). As mentioned in [15], “ Watch out for frequency 
warping if the filter bandwidth is larger than about 20% of the simulation bandwidth (1/sampling 
period);; the passband frequencies will be correct, but out-of-band characteristics may be warped.” 
Thus, to avoid frequency warping, the following condition has to be met: 

&r<(Q2)-f,\ 

where BT = Bandwidth-Time product (defined in Appendix A); and 
f s = sampling frequency. 

For simulations involving BT=1, f s =16 was used and for BT=2,3 values of f s =32, 48 were used to 
avoid frequency warping. But what causes frequency warping? As described in [10], frequency 
warping is caused by the nonlinear relationship between corresponding digital and analog frequencies. 
Therefore if the sampling frequency is not high enough the high frequencies for the digital filter will 
be warped compared to its corresponding analog filter. 
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APPENDIX D 


PROGRAM LISTING: 

“SPWSOBRE.M” 

This appendix contains the listing of the program SPWSOBRE.M written in Matlab that was used 
to produce the scatter plots output from the Non-Constant Envelope simulations. It calculates the 
Average Symbol Variance of each of the 8 decision regions for different values of BT. 


Program Listing: SPWSOBRE.M 

% Program : SPWSOBRE.M 
% Written by : Ruben Caballero 

% Developed using MATLAB V4.2B and Signal Processing Toolbox V3.0b 
% Copyright 1 996 - Ruben Caballero 
% New Mexico State University 

% 

% This function loads the file from the non-constant envelope simulations performed 
% in SPW, and locates the decision region where each symbol belongs. Afterwards, 
% the function will calculate the mean and variance of each decision region using 
% the symbols. For each simulation, 625 symbols were used to produce the output 
% scatter plots and calculate the average symbol variance. 


o /o *«****«*******»«* INPUTS; CHANGE THESE VARIABLES*************** 

clear; 

tic; 

load sambutl7.asc; % file containing the symbols 


jj=i; 

for ii=38: length(sambut 17); % 38 specifies the filter delay: change this value depending ... 

f° r iipl ; 2 % ... on the simulation and filter 

A(ijjjj)=sambutl7(iijjj); % Change also the name of the ASCII file 
end; 

jj=jj+i; 

end; 


nsymbol=48; 


% sampling frequency to be changed depending on the ... 
% the simulation 
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0 /q ************* ***********pj^QQp^jy| STARTS************************ 


[nr nc]=size(A); 
for q=l:nr 

x(q) = A(q, 1 ) * (cos( A(q,2))+sin( A(q,2))* i) ; 
end; 

[nr nc]=size(x); 
Ksymbols=floor(nc/nsymbol); 
x=x( 1 :Ksymbols*nsymbol); 


z=[-l:0.01:l]; % Create the axis for the Scatter Plot 

[rz nz]=size(z); 

y=sqrt(ones( 1 ,nz)-z. * z); 

plotfey/b-'^-y/b-’); 

grid 

axis('square'); 
hold on; 

avg=sum(abs(x)); 

avg=avg/nc; 


for i=l:Ksymbols 
m=(nsymbol)*i; 
symbols(i)=x(m)/avg; 
end 

% Determine the location of the symbols. 

RxSig=angle(symbols); 

decO=l; 
decl=l; 
dec2=l; 
dec3=l; 
dec4=l; 
dec5=l; 
dec6= 1 ; 
dec7=l ; 

for i= 1 :length(RxSig) % Decision Region 0 

if ((RxSig(i)>=0) & (RxSig(i) <=pi/4)) 
decisionO(decO)=symbols(i); 
decO=decO+l ; 

end; 
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% Decision Region 1 


if ((RxSig(i)>=pi/4) & (RxSig(i) <=pi/2)) 
decision 1 (dec 1 )=symbols(i); 
decl=decl + l; 

end; 

if ((RxSig(i)>pi/2) & (RxSig(i) <=6*pi/8)) % Decision Region 2 

decision2(dec2)=symbols(i); 
dec2=dec2+l; 

end; 

if ((RxSig(i)>=6*pi/8) & (RxSig(i) <=pi)) % Decision Region 3 

decision3(dec3)=syrnbols(i); 
dec3=dec3+l; 

end; 


if ((RxSig(i)>-pi) & (RxSig(i) <= -6* pi/8)) % Decision Region 4 

decision4(dec4)=symbols(i); 
dec4=dec4+l; 

end; 

if ((RxSig(i)>-6*pi/8) & (RxSig(i) <=-pi/2)) % Decision Region 5 

decision5(dec5)=symbols(i); 
dec5=dec5+l; 

end; 

if ((RxSig(i)> -pi/2) & (RxSig(i) <=-pi/4)) % Decision Region 6 

decision6(dec6)=symbols(i); 
dec6=dec6+l; 

end; 

if ((RxSig(i)>-pi/4) & (RxSig(i) <0)) % Decision Region 7 

decision7(dec7)=symbols(i); 
dec7=dec7+l; 

end; 

end; % end FOR 
pio^symbols/go'); 
axis([-1.25 +1.25 -1.25 +1.25]); 

title('Output of SSPA (saturation level), 5th Order Butterworth Filter (BT=3.8)’); 

xlabel(’rear); 

ylabel('imag'); 

% Calculation of the Mean of the Symbols in the Decision Regions 

meanrad(l)=(mean(angle(decisionO))); % mean in rads 
meandeg( 1 )=(meanrad( 1 )* 1 80/pi); % mean in degrees 

meanrad(2)=(mean(angle(decision 1 ))); % mean in rads 

meandeg(2)=(meanrad(2)* 1 80/pi); % mean in degrees 
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meanrad(3)=(mean(angle(decision2))); 
meandeg(3)=(meanrad(3)* 180/pi); 

% mean in rads 
% mean in degrees 

meanrad(4)=(mean(angle(decision3))); 
meandeg(4)=(meanrad(4)* 1 80/pi); 

% mean in rads 
% mean in degrees 

meanrad(5)=(mean(angle(decision4))); 
meandeg(5)=(nieanrad(5)* 1 80/pi); 

% mean in rads 
% mean in degrees 

meanrad(6)=(mean(angle(decision5))); 
meandeg(6)=(meanrad(6)* 1 80/pi); 

% mean in rads 
% mean in degrees 

meanrad(7)=(mean(angle(decision6))); 
meandeg(7)=(meanrad(7)* 1 80/pi); 

% mean in rads 
% mean in degrees 

meanrad(8)=(mean(angle(decision7))); 
meandeg(8)=(meanrad(8)* 1 80/pi); 

% mean in rads 
% mean in degrees 

% Calculation of the Variance of the Symbols in the Decision Regions 

varrad( 1 )=(std(angle(decision0))). A 2; 
vardeg( 1 )=(varrad( 1 )* 1 80/pi); 

% variance in rads 
% variance in degrees 

varrad(2)=(std(angle(decision 1 ))). A 2; 
vardeg(2)=(varrad(2)* 1 80/pi); 

% variance in rads 
% variance in degrees 

varrad(3) : =(std(angle(decision2))). A 2; 
vardeg(3)=(varrad(3)* 1 80/pi); 

% variance in rads 
% variance in degrees 

varrad(4)=(std(angle(decision3))). A 2; 
vardeg(4)=(varrad(4)* 1 80/pi); 

% variance in rads 
% variance in degrees 

varrad(5)=(std(angle(decision4))) . A 2 ; 
vardeg(5)=( varr ad(5)* 1 80/pi); 

% variance in rads 
% variance in degrees 

varrad(6)=(std(angle(decision5))). A 2; 
vardeg(6)=(varrad(6)* 1 80/pi); 

% variance in rads 
% variance in degrees 

varrad(7)=(std(angle(decision6))). A 2; 
vardeg(7)=(varrad(7)* 1 80/pi); 

% variance in rads 
% variance in degrees 

varrad(8)=(std(angle(decision7))). A 2; 
vardeg(8)=(varrad(8)* 1 80/pi); 

% variance in rads 
% variance in degrees 


save meandeg; save vardeg; 
time=(toc/60); % in minutes 
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APPENDIX E 


SCATTER PLOTS: 

5 th Order BUTTERWORTH FILTER 

This appendix contains the Scatter Plots of the output of the SSPA when a 5th Order Butterworth 
Spectrum Filter is used after the 8PSK Modulator. The Scatter Plots were produced for BT=1 to 4 
using the SPWSOBRE.M function. Also, the listing of the program BUTTPLOT.M written in Matlab 
which was used to produce the Average Symbol Variance plots (see Figure 3.43 and 3.44 in this 
report) is included. All the values of mean and variance in this program were calculated in the 
program SPWSOBRE.M. 

E.l Scatter Plots 
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Figure E.2 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: 5th Order Butterworth BT=1.05 
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- Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: 5th Order Butterworth BT=3.8 
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E.2 Program Listing: BUTTPLOT.M 


% Program : BUTTPLOT.M 
% Written by : Ruben Caballero 

% Developed using MATLAB V4.2B and Signal Processing Toolbox V3.0b 
% Copyright 1 996 - Ruben Caballero 
% New Mexico State University 
% 

% This function uses the mean and variance calculated by the function SPWSOBRE.M to 
% calculate the Average Symbol Variance. The result is a plot of the Average Symbol Variance 
% versus BT for the 5th Order Butterworth filter as the spectrum shaping filter. 

clear; 


% ***************** Butterworth Filters ******************************** 


BT=[1 1.05 1.1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3,2 3.4 3.6 3.8 4.0]; % set BT values 

%* *** * ****** * **** * **** * *** * ********** 4c **** * ******************* *** * * 

%* * * * * BT= 1 fs =16 ************************************************* 

%* **************************************************************** * 

% BUTTERWORTH: BT=1 (using delay =20, ascii file: sampbutl.asc) 

meandeg(l,:)= [1.7020e+001 6.2250e+001 1.0680e+002 1.5258e+002 -1.6298e+002 
-1.1786e+002 -7.2650e+001 -2.7753e+001]; 

vardeg(l,:)= [7.1477e-002 6.3547e-002 5.2001e-002 6.6354e-002 5.4146e-002 5.4370e-002 
6.6427e-002 7.4719e-002]; 
avgvardegl = sum(vardeg(l,:))/8; 


<^****************************************************************** 

%***♦* BT=1.05 fs =16 ********************************************** 

%* *** * ******** * ********* * ** * * * * * * * * * * * * * * * * * * * ** * * * * * * * * * * 4c 4c * * * * * * * 

% BUTTERWORTH: BT=1.05 (using delay =13, ascii file: sampbutl3.asc) 

meandeg(13,:) = [1.7092e+001 6.2220e+001 1.0703e+002 1.5225e+002 -1.6294e+002 
-1.1 794e+002 -7.2843e+001 -2.7899e+001]; 

vardeg(13,:) = [1.6673e-002 1.3528e-002 1.2779e-002 1.4959e-002 1.2526e-002 1.2054e-002 
1 ,3539e-002 1.5471e-002]; 
avgvardeg 1 3 = sum(vardeg(13,:))/8; 
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%* ************** * ******** * * * * * * * * * * ]|c * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 

%***** BT=1.1 fs =32 *********************************************** 

%* ***********************************************************^*^ )>jJ( 

% BUTTERWORTH: BT=1.1 (using delay =24, ascii file: sampbut8.asc) 

meandeg(8,:) = [1.7103e+001 6.2082e+001 1.0716e+002 1.5204e+002 -1.6296e+002 
-1.1789e+002 -7.2971e+001 -2.7915e+001]; 

vardeg(8,:) = [1.8056e-003 1.8445e-003 1.9903e-003 2.0978e-003 2.0085e-003 2.1336e-003 
1.7640e-003 1.8487e-003]; 
avgvardeg8 = sum(vardeg(8,:))/8; 

0 ^* **************************************************************** 

%****♦ BT=1.2 fs =16 *********************************************** 

O/^* **************************************************************** 

% BUTTERWORTH: BT=1.2 (using delay =18, ascii file: sampbut4.asc) 

meandeg(4,:)= [1.7096e+001 6.2134e+001 1.0702e+002 1.5222e+002 -1.6285e+002 
-1.1 789e+002 -7.2802e+001 -2.7860e+001]; 

vardeg(4,:) = [2.7145e-003 2.6749e-003 2.5493e-003 2.6131e-003 2.0687e-003 2.9801e-003 
2.5173e-003 2.3366e-003]; 
avgvardeg4 = sum(vardeg(4,:))/8; 

%* **************************** ************************************ 

%***** 3^=1 4 fs =32 *********************************************** 
%*********************♦♦*************************♦**************** 

% BUTTERWORTH: BT=1.4 (using delay =19, ascii file: sampbutS.asc) 

meandeg(5,:)= [1.7092e+001 6.2088e+001 1.0717e+002 1.5205e+002 -1.6299e+002 
-1.1789e+002 -7.2985e+001 -2.7915e+00I]; 

vardeg(5,:) = [3.6223e-003 3.0185e-003 3.5466e-003 3.6122e-003 3.3108e-003 3.5702e-003 
3.4100e-003 2.9896e-003]; 
avgvardegS = sum(vardeg(5,:))/8; 

o^******************************* *************************** ******* 

%***** BT=1.6 fs =32 *********************************************** 

o^****************************** *********************************** 


% BUTTERWORTH: BT=1,6 (delay =27, ascii file: sampbut6.asc) 

meandeg(6,:)==[1.7105e+001 6.2079e+001 1.0712e+002 1.5207e+002 -1.6290e+002 
-1.1789e+002 -7.2922e+001 -2.790 le+001]; 

vardeg(6,:) = [5.7935e-004 4.9844e-004 4.7372e-004 5.3648e-004 4.6827e-004 5.2152e-004 
4.4754e-004 4.9717e-004]; 
avgvardeg6 = sum(vardeg(6,:))/8; 
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%** ******************************************************** ******* 

°/ q * * * * * 3T=1.8 fs=48 *********************************************** 

o / Q * *************************************************************** * 

% BUTTERWORTH: BT=1.8 (delay =35, ascii file: sampbut7.asc) 

meandeg(7,:)= [1.7103e+001 6.2156e+001 1.0708e+002 1. 52 17e+002 -1.629 le+002 
-1.1791e+002 -7.287 le+001 -2.7893e+001]; 

vardeg(7,:) = [3.0054e-003 2.4409e-003 2.3048e-003 2.7448e-003 2.2233e-003 2.1648e-003 
2.4702e-003 2.8097e-003]; 
avgvardeg7 = sum(vardeg(7,:))/8; 

%* *************************************************************** * 
o^* * * * * 3T— 2 fs=48 ************************************************* 

%* *************************************************************** * 

%BUTTERWORTH: BT=2 (using delay = 32, ascii file: sampbut2.asc) 

meandeg(2,:) = [1.7100e+001 6.211 le+001 1.0710e+002 1.521 le+002 -1.6290e+002 - 
1.1 790e+002 -7.2897e+001 -2.7899e+001]; 

vardeg(2,:) = [1.2757e-004 1.0112e-004 9.9469e-005 1.1625e-004 9.3600e-005 8.5729e-005 
1.0797e-004 1 .221 8e-004]; 
avgvardeg2 = sum(vardeg(2,:))/8; 


o / 0 * *************************************************************** * 

^***** gy =2 2 fs=48 ************************************************ 

o / 0 * *************************************************************** * 

%BUTTERWORTH: BT=2.2 (using delay = 29, ascii file: sampbut9.asc) 

meandeg(9,:) = [1.7098e+001 6.2134e+001 1.0707e+002 1.5215e+002 -1.6289e+002 
-1 .179 le+002 -7.2870e+001 -2.7895e+001]; 

vardeg(9,:)= [1.101 9e-003 9.3545e-004 8,7686e-004 1.0168e-003 8.5896e-004 9.3158e-004 
8.5970e-004 9.7456e-004]; 
avgvardeg9 = sum(vardeg(9,:))/8; 


o^***** gj = 2 4 fs=48 *********************************************** 


%BUTTERWORTH: BT=2,4 (using delay = 27, ascii file: sambutlO.asc) 

meandeg(10,:) = [1.7080e+001 6.2064e+001 1.0709e+002 1.5206e+002 -1.6292e+002 
-1.1 792e+002 -7.2927e+00 1 -2.7922e+00 1 ]; 

vardeg(10,:) = [2.5908e-004 2.1067e-004 2.0091e-004 2.3544e-004 1.9016e-004 1.8373e-004 
2.0921e-004 2.4034e-004]; 
avgvardeglO = sum(vardeg(10,:))/8; 
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%* ***************************************** * ** * *** * * * * * * * * * * * * * * * * 

%***** 37=2 6 fs=48 ************************************************ 

% ********************* * * * * * * * * * * * * 4c * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 

%BUTTERWORTH: BT=2.6 (using delay =34, ascii file: sambutl l.asc) 

meandeg(ll,:) = [1.7076e+001 6.2055e+001 1.0709e+002 1.5205e+002 -1.6293e+002 
-1.1792e+002 -7.2940e+001 -2.7928e+001]; 

vardeg(l 1,:) = [3.9346e-004 3.3146e-004 3.0956e-004 3.6135e-004 3.0312e-004 3.1763e-004 
3.0654e-004 3.5248e-004]; 
avgvardegl 1 = sum(vardeg(l l,:))/8; 

% ******************** * ***** * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 

%*♦*** 37=2.8 fs=48 *********************************************** 

o^********************************************************^*^^^ 


%BUTTERWORTH: BT-2.8 (using delay =32, ascii file: sambut!2.asc) 

meandeg(12,:) = [1.7078e+001 6.2083e+001 1.0708e+002 1.5209e+002 -1.6292e+002 
-1.1 792e+002 -7.2917e+001 -2.792 le+001]; 

vardeg(12,:) = [1.4141e-005 1.2143e-005 1.2144e-005 1.3224e-005 l.I092e-005 1.2642e-005 
1.1685e-005 1 ,2637e-005]; 
avgvardegl2 = sum(vardeg(12,:))/8; 

%* ******* **************** * * * * * * * * * * * * * * ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ # ^ # + + ^ # + + * * * 

%***** BT=3 fs=48 ************************************************* 

%* ********************** 4c ***** * ** * * * * * * * * * * * * * * * * * 4c * * * * 41 * * * * * * * * * * 

%BUTTER W ORTH : BT=3 (using delay =30, ascii file: sampbut3.asc) 

meandeg(3,:) = [1.7077e+001 6.2088e+001 1.0707e+002 1.5209e+002 -1.6292e+002 
-1.1792e+002 -7.2913e+001 -2.792 le+001]; 

vardeg(3,:) = [9.9154e-005 8.3102e-005 7.9266e-005 9.1259e-005 7.4914e-005 8.1407e-005 
7.9067e-005 8.9380e-005]; 
avgvardeg3 = sum(vardeg(3,:))/8; 

o^* ************** ************ ************************************** 

%***** 37=3.2 fs=48 ************************************************ 
0 ^*************** ******************************************** + # 


%BUTTERWORTH: BT-3.2 (using delay = 28, ascii file: sambut!4.asc) 

meandeg(14,:) = [1.7072e+001 6.2073e+001 1.0707e+002 1.5207e+002 -1.6293e+002 
-1.1793e+002 -7.2927e+001 -2.7928e+001j; 

vardeg(14,:)= [5.8964e-007 5.0913e-007 9.2412e-007 6.3686e-007 6.4852e-007 4.5768e-007 
1 0022e-006 5.8489e-007]; 
avgvardegl4 = sum(vardeg(14,:))/8; 
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%********* ****** ************************************************** 

% ***** BT=3.4 fs=48 ************************************************ 

0/ 0 * ******** ^****^ic*^*****t ******************************* ********** 

%BUTTERWORTH: BT=3.4 (using delay = 34, ascii file: sambutlS.asc) 

meandeg(15,:) = [1.7069e+001 6.2062e+001 1.0707e+002 1.5206e+002 -1.6293e+002 
-l.I793e+002 -7.2936e+001 -2.7932e+001]; 

vardeg(15,:)=[4.8152e-005 4.0562e-005 3.7986e-005 4.4103e-005 3.8410e-005 3.8316e-005 
3.7664e-005 4.32 10e-005]- 
avgvardegl5 = sum(vardeg(15,:))/8; 

0 ,/ ************************************************ ***************** 

o/ o ***** 37=3 5 fs=48 *********************************************** 
o/ 0 ** ******** * ** ** * **************** ** ** **** * ************************ 


%BUTTERWORTH: BT=3.6 (using delay = 32, ascii file: sambutl 6 .asc) 

meandeg(16,:)=[1.7068e+001 6.2064e+001 1.0707e+002 1.5206e+002 -1.6293e+002 
-1.1793e+002 -7.2934e+001 -2.7933e+001]; 

vardeg(16,:) = [1.6145e-005 1.3681e-005 1.2928e-005 1.4785e-005 1.3520e-005 1.2891e-005 
1.2778e-005 1.443 le-005]; 
avgvardegl 6 = sum(vardeg(16,:))/8; 


< y o ***** 37=3 g fs=48 *********************************************** 


% BUTTERWORTH: BT=3.8 (using delay =37, ascii file: sambutl7.asc) 

meandeg(17,:)= [1.7068e+001 6.2071e+001 1.0707e+002 1.5207e+002 -1 .6293e+002 
-1.1793e+002 -7,2929e+001 -2.7932e+001]; 

vardeg(17,:)= [6.081 8e-006 5.0583e-006 5.1543e-006 5.6452e-006 4.5621e-006 4.9250e-006 
5.5862e-006 5.5223e-006]; 
avgvardegl7 = sum(vardeg(17,:))/8; 




% BUTTERWORTH: BT=4.0 (using delay =35, ascii file: sambutl 8 .asc) 

meandeg(18,:)= [I.7067e+001 6.2068e+001 1.0707e+002 1.5207e+002 -1.6293e+002 
-1.1793e+002-7.2932e+001 -2.7933e+001]; 

vardeg( 1 8,:) = [6.3 1 59e-007 5.2962e-007 7.6612e-007 6.0844e-007 7.0038e-007 5.1824e-007 
1.1054e-006 5.8058e-007]; 
avgvardegI 8 = sum(vardeg(18,:))/8; 


149 



%**************************************** 

o/ o ** **************** PLOTS*** ************ 
0 / 0 * *************************************** 


% Average Variance for the Symbols Plot 

figure(l); 

whitebg; 

AMPMavgvar = [avgvardegl avgvardegl3 avgvardeg8 avgvardeg4 avgvardeg5 avgvardeg6 
avgvardeg7 avgvardeg2 avgvardeg9 avgvardeglO avgvardegl 1 avgvardegl 2 
avgvardeg3 avgvardegl 4 avgvardegl 5 avgvardegl 6 avgvardegl 7 avgvardegl 8]; 

% interpolation using spline fit 

interBT= 1:0.0 1:4.0; 

interpol=spline(BT,AMPMavgvar,interBT); 

plot(BT,AMPMavgvar,'ok',interBT,interpol,'k'); 

title(’5th Order Butterworth Filter Symbol Average Variance vs BT); 
xlabel('Bandwidth-Time Product (BT)'); 
ylabel(' Variance in Degrees'); 
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APPENDIX F 


SCATTER PLOTS: 

3 rd Order BESSEL FILTER 

This appendix contains the Scatter Plots of the output of the SSPA when a 3rd Order Bessel 
Spectrum Filter is used after the 8PSK Modulator. The Scatter Plots were produced for BT=1 to 4 
using the SPWSOBRE.M function. Also, the listing of the program BESSPLOT.M written in Matlab 
which was used to produce the Average Symbol Variance plots (see Figure 3.45 and 3.46 in this 
report) is included. All the values of mean and variance in this program were calculated in the 
program SPWSOBRE.M. 

F.l Scatter Plots 
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re F.7 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: 3rd Order Bessel BT= 1 .8 
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Figure F.8 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: 3rd Order Bessel BT=2 
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Output of SSPA (saturation level), 3rd Order Bessel Filter (BT=2.4) 
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Figure F. 1 0 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: 3rd Order Bessel BT=2.4 
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Output of SSPA (saturation level), 3rd Order Bessel Filter (BT=2.8) 
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Figure F. 12 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: 3rd Order Bessel BT=2.8 
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Figure F.15 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: 3rd Order Bessel BT=3. 
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Output of SSPA (saturation level), 3rd Order Bessel Filter (BT=4) 



real 


Figure F. 1 8 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: 3rd Order Bessel BT=4 
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F.2 Program Listing: BESSPLOT.M 

% Program : BESSPLOT.M 
% Written by : Ruben Caballero 

% Developed using MATLAB V4.2B and Signal Processing Toolbox V3.0b 
% Copyright 1996 - Rub£n Caballero 
% New Mexico State University 
% 

% This function uses the mean and variance calculated by the function SPWSOBRE.M to 
% calculate the Average Symbol Variance. The result is a plot of the Average Symbol Variance 
% versus BT for the 3rd Order Bessel filter as the spectrum shaping filter. 

clear; 

<y o * **** * *********** ** ***************** * ***** * * * **** * * * ********* * *** * 

%******* ******* *** Bessel Filters ************************************* 

0/ o ****************************************************************** 


BT=[1 1.05 1.1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4.0]; 


o / Q * ***************************************************************** * 

%***♦* 3 T=] fs=16 *************************************************** 

o/ 0 * ****************************************************************** 


% BESSEL : BT=1 (using delays 10, ascii file: sampbesl.asc) 

meandeg(l,:) = [1.7100e+001 6.2067e+001 1.0712e+002 1.5205e+002 -1.6290e+002 
-1.1789e+002 -7.2926e+001 -2.7905e+001]; 

vardeg(l,:) = [1.022 le-003 8.5808e-004 8.0633e-004 9.4189e-004 7.8621e-004 8.2393e-004 
8.1 1 16e-004 9.1752e-004]; 
avgvardegl = sum(vardeg(l,:))/8; 

oy Q * ***************************************************************** * 

&/ q * * * * * 05 fs=32 ************************************************ 

<y o ***** ******************************************************** ****** 


% BESSEL : BT=1.05 (using delay =20, ascii file: sampbes!3.asc) 

meandeg(13,:) = [1.7098e+001 6.2105e+001 1.0709e+002 1.521 le+002 -1.6290e+002 
-1.1790e+002 -7.2897e+001 -2.790 le+001]; 

vardeg(13,:)= [4.0437e-005 3.3809e-005 3.1692e-005 3.7052e-005 3.0734e-005 3.2129e-005 
3.21 10e-005 3.6434e-005]; 
avgvardegl 3 = sum(vardeg(13,:))/8; 
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0 / 0 * ********** * ******************* * ***** * * * * * * + ^ + + ^ ^ + * ^ ^ ^ ^ ^ ^ # ^ ^ # ^ ^ ^ 

%**♦** BT=1.1 fs=32 ************************************************* 

%* ***** * * *** * ******** * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 

% BESSEL : BT=1. 1 (using delay = 19, ascii file: sampbes8.asc) 

meandeg(8,:) = [1.7098e+001 6.2101e+001 1.0710e+002 1.5210e+002 -1.6290e+002 
-1.1790e+002 -7.2900e+001 -2.7902e+001]; 

vardeg(8,:)=[9.9913e-006 8.3489e-006 7.8251e-006 9.1203e-006 7.5720e-006 7.9193e-006 
7.9803e-006 9.041 le-006]; 
avgvardeg8 = sum(vardeg(8,:))/8; 

%*********** ******************************** ************************ 

%****♦ BT=1.2 fs=32 ************************************************* 
%***^************* *********************************** ^^^^^^^ ** 

% BESSEL : BT=1.2 (using delay =17, ascii file: sampbes4.asc) 

meandeg(4,:)= [1.7099e+001 6.2083e+001 1.071 le+002 1.5208e+002 -1.6290e+002 
-1.1790e+002 -7.2914e+001 -2.7903e+001]; 

vardeg(4,:) = [2.4972e-004 2.0951e-004 1.9659e-004 2.3032e-004 1.9168e-004 2.0041e-004 
1.9809e-004 2.2399e-004]; 
avgvardeg4 = sum(vardeg(4,:))/8; 




% BESSEL : BT-1.4 (using delay = 24, ascii file: sampbes5.asc) 

meandeg(5,:) = [1.7101e+001 6.2102e+001 1.0710e+002 1.5210e+002 -1.6290e+002 
-1.1790e+002 -7.2899e+001 -2.7899e+001]; 

vardeg(5,:)= [1.3749e-006 1.1572e-006 1.1217e-006 1.2701e-006 1.1957e-006 1.2262e-006 
1.1256e-006 1.2264e-006]; 
avgvardeg5 = sum(vardeg(5,:))/8; 


%******* M ************«** , '********«******************««**********»*« 

%***** BT=1 .6 fs=32 ************************************************* 
%******************************************************************* 

% BESSEL : BT=1 .6 (using delay = 22, ascii file: sampbes6.asc) 

meandeg(6,:) = [ 1 .71 06e+00 1 6.2104e+001 1.0711e+002 1.5210e+002 -1.6289e+002 
-1.1789e+002 -7.2896e+001 -2.7894e+001]; 

vardeg(6,:) = [4.5746e-006 3.8329e-006 3.5877e-006 4.2036e-006 3.3338e-006 3.5724e-006 
3.6400e-006 4.1 172e-006]; 
avgvardeg6 = sum(vardeg(6,:))/8; 
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%* ************************************ * *** * * * * * * *** * * * * * * * * * * * * * * * * * 
%****** BT=1.8 fs=48 ************************************************ 

o^** ********* **********41 ********************************************* 


% BESSEL : BT=1.8 (using delay = 29, ascii file: sampbes7.asc) 

meandeg(7,:) = [1.7109e+001 6.2107e+001 1.071 le+002 1.521 le+002 -1.6289e+002 
-1.1789e+002 -7.2892e+001 -2.789 le+001]; 

vardeg(7,:) = [1.6004e-006 1.3463e-006 1.2807e-006 1.4545e-006 1.1838e-006 1.2862e-006 
1.3147e-006 1 ,4452e-006]; 
avgvardeg7 = sum(vardeg(7,:))/8; 


%* ****************************************************************** 

%****** BT=2 fs=32 ************************************************** 

O^j* ****************************************************************** 

% BESSEL : BT=2 (using delay = 17, ascii file: sampbes2.asc) 

meandeg(2,:) = [1.7109e+001 6.2109e+001 1.071 le+002 1.521 le+002 -1.6289e+002 
-1.1 789e+002 -7.289 le+001 -2.789 le+001]; 

vardeg(2,:) = [4.4990e-008 3.8308e-008 5.1507e-008 5.3124e-008 5.0042e-008 5.9081e-008 
6.3786e-008 4.0034e-008]; 
avgvardeg2 = sum(vardeg(2,:))/8; 


%* *************************************************** * * * * * * * * * * * * * * * 
%♦**** BT=2.2 fs=48 ************************************************* 
%* ************************************************** * * * * * * * * * * * * * * * * 


% BESSEL : BT=2.2 (using delay = 31, ascii file: sampbes9.asc) 

meandeg(9,:)=[1.7108e+001 6.2106e+001 1.071 le+002 1.521 le+002 -1.6289e+002 -1.1 789e+002 
-7.2893e+001 -2.7893e+001]; 

vardeg(9,:) = [1.0595e-006 8.9530e-007 8.6564e-007 9.7955e-007 9.1845e-007 9.8303e-007 
8.653 le-007 9.4939e-007]; 
avgvardeg9 = sum(vardeg(9,:))/8; 

%* ********* ********************************************************* 

%****♦ BT— 2.4 fs=48 ************************************************* 

%* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 4c * * * ale * * * * * 


% BESSEL : BT=2.4 (using delay = 29, ascii file: sambeslO.asc) 

meandeg(10,:) = [1.7105e+001 6.2104e+001 1.0711e+002 1.5210e+002 -1.6290e+002 
-1.1790e+002 -7.2896e+001 -2.7895e+001]; 

vardeg(10,:) = [5.5214e-007 4.6563e-007 4.7792e-007 5.3915e-007 5.5696e-007 5.8688e-007 
4.5275e-007 4.8757e-007]; 
avgvardeglO = sum(vardeg(10,:))/8; 
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% BESSEL : BT=2.6 (using delay = 27, ascii file: sambes 1 1 .asc) 

meandeg( 1 1 ,:) = [ 1 .7 1 0 1 e+00 1 6.2100e+001 1.0710e+002 1.5210e+002 -1.6290e+002 
-1.1 790e+002 -7.2899e+001 -2.7899e+001]; 

vardeg(l 1,:) = [4.0654e-007 3.4041e-007 3.8516e-007 4.1756e-007 4.7282e-007 4.7462e-007 
3.4439e-007 3.5223e-007]; 
avgvardegl 1 = sum(vardeg(l l,:))/8; 

%****************************************************************,», 

%***** BT=2.8 fs=48 ************************************************* 

"/a***************************************************************,**** 


% BESSEL : BT=2.8 (using delay = 25, ascii file: sambesl2.asc) 

meandeg(12,:) = [1.7096e+00I 6.2095e+001 1.0710e-002 1.5209e+002 -1.6290e+002 
-1.1 790e+002 -7.2904e+001 -2.7904e+001]; 

vardeg(12,:) = [4.5417e-007 3.8038e-007 4.5201e-007 4.5872e-007 5.5238e-007 5.1875e-007 
4.0865e-007 3.9254e-007]; 
avgvardegl 2 = sum(vardeg(12,:))/8; 


% BESSEL : BT=3 (using delay = 40, ascii file: sampbesS.asc) 

meandeg(3,:)= [1.7092e+001 6.2092e+001 1.0709e+002 1.5209e+002 -1.629 le+002 
-1.1791e+002 -7.2908e+001 -2.7908e+001]; 

vardeg(3,:)= [5.829 le-008 3.9094e-008 1.2009e-007 1.2832e-007 1.0620e-007 1.1 197e-007 
7. 1450e-008 3.9674e-008]; 
avgvardeg3 = sum(vardeg(3,:))/8; 


% BESSEL : BT=3.2 (using delay = 28, ascii file: sambesM.asc) 

meandeg(14,:)= [1.7088e+001 6.2088e+001 1.0709e+002 1.5209e+002 -1.629 le+002 
-1.1 79 le+002 -7.291 1 e+00 1 -2.7912e+001]; 

vardeg(14,:)= [5.4259e-008 3.5312e-008 1.4483e-007 1.0340e-007 I.1655e-007 8.1064e-008 
1.3418e-007 4.1300e-008]; 
avgvardeg!4 = sum(vardeg(14,:))/8; 
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%* ****************************************************************** 
%*♦*** BT=3.4 fs=48 ************************************************* 
<y 0 * ****************************************************************** 


% BESSEL : BT=3.4 (using delay = 28, ascii file: sambesl5.asc) 

meandeg( 1 5 , :) = [ 1 . 7084e+00 1 6.2084e+001 1.0708e+002 1.5208e+002 -1.6292e+002 
-1.1792e+002 -7.2915e+001 -2.7916e+001]; 

vardeg(15,:) = [4.0005e-008 2.5838e-008 1.4176e-007 7.8527e-008 1.1756e-007 6.8254e-008 
1.4721e-007 2.9590e-008]; 
avgvardegl5 = sum(vardeg(15,:))/8; 

%************************************************* ********** ******** 

%****♦ BT=3 6 fs=48 ************************************************* 

%* * ******* * ********* ** ** * * * **** * * ****** ****** * ** * * « * * * ***** * ** * * * * * * 


% BESSEL : BT=3.6 (using delay = 28, ascii file: sambesl6.asc) 

meandeg(16 ? :) = [ 1.708 le+001 6.2081e+001 L0708e+002 L5208e+002 -1.6292e+002 
-1.1792e+002 -7.2919e+001 -2.7919e+001]; 

vardeg(16,:) = [3.1934e-008 2.0872e-008 1.3377e-007 6.1400e-008 1.0893e-007 5.5746e-008 
1.5336e-007 2.3636e-008]; 
avgvardegl6 = sum(vardeg(16,:))/8; 


o / o * ****************************************************************** 

<^/***** BT=3 S f s =48 ************************************************* 
0 / 0 * ****************************************************************** 


% BESSEL : BT=3.8 (using delay = 28, sambesl7.a$c) 

meandeg(17,:) = [L7078e+001 6.2078e-K)01 1.0708e+002 L5208e+002 -1.6292e+002 
-1.1792e+002 -7.2921e+001 -2.7922e+001]; 

vardeg(17,:) = [2.5426e-008 L6786e-008 1.2010e-007 4.7555e-008 9.4054e-008 4.3259e-008 
1.5054e-007 1.9058e-008]; 
avgvardegl7 = sum(vardeg(17,:))/8; 

%* ****************************************************************** 
o^***** fs=4g *************************************************** 

o / o * ****************************************************************** 


% BESSEL : BT=4 (using delay = 28, sambesl8.asc) 

meandeg(18,:)= [1.7076e+001 6.2075e+001 1.0708e+002 L5208e+002 -1.6292e+002 
-1.1792e+002 -7.2924e+001 -2.7924e+001]; 

vardeg(18,:)= [1.9867e-008 L3236e-008 1.0380e-007 3.6088e-008 7.8051e-008 3.2460e-008 
1 .4042e-007 1.5133e-008]; 
avgvardeg!8 = sum(vardeg(18,:))/8; 
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%************** pLQyg*************** 


% Average Variance for the Symbols 

figure(l); 

whitebg; 

AMPMavgvar = [avgvardegl avgvardegl3 avgvardeg8 avgvardeg4 avgvardegS avgvardeg6 

avgvardeg7 avgvardeg2 avgvardeg9 avgvardeglO avgvardegl 1 avgvardegl 2 
avgvardeg3 avgvardegl 4 avgvardegl 5 avgvardegl 6 avgvardegl 7 avgvardegl 8]; 

% interpolation using spline fit 

interBTl = l :0. 1:4.0; 

interpol 1 =interp 1 (BT,AMPMavgvar,interBTl 'spline'); 
plot(BT,AMPMavgvar,'ok’,interBT 1 , (interpol 1 )/k'); 

plotCBT^MPMavgvar/ok’); 

title('3rd Order Bessel Filter Symbol Average Variance vs BP); 
xlabel(’Bandwidth-Time Product (BT)'); 
ylabel(' Variance in Degrees 1 ); 
grid; 
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APPENDIX G 


SCATTER PLOTS: 

SRRC FILTER 

This appendix contains the Scatter Plots of the output of the SSPA when a Square Root Raised 
Cosine (SRRC) Filter is used after the 8PSK Modulator. The Scatter Plots were produced for a=0. 1 to 
1 using the SPWSOBRE.M function. Also, the listing of the program SRRCPLOT.M written in 
Matlab which was used to produce the Average Symbol Variance plot (see Figure 3.47 in this report) 
is included. All the values of mean and variance in this program were calculated in the program 
SPWSOBRE.M. 

G.l Scatter Plots 
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•e G.l - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: SRRC a = 0.1 
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re G.3 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: SRRC a=0.2 
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- Output of SSPA (Saturation Level) 
Spectrum Shaping Filter: SRRC a=0.25 
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e G.8 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: SRRC a=0.45 
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igure G.9 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: SRRC a=0.5 
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0 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: SRRC a=0.55 
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re G. 1 1 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: SRRC a=0.6 
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Figure G. 12 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: SRRC a=0.65 
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imag 


Output of SSPA (saturation level), SRRC Filter (alpha=0.7) 



Figure G. 13- Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: SRRC a=0.7 
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- Output of SSPA (Saturation Level) 
Spectrum Shaping Filter: SRRC a=0.75 
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Figure G. 1 8 - Output of SSPA (Saturation Level) 

Spectrum Shaping Filter: SRRC a=0.95 
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G.2 Program Listing: SRRCPLOT.M 

% Program : SRRCPLOT.M 
% Written by : Ruben Caballero 

% Developed using MATLAB V4.2B and Signal Processing Toolbox V3.0b 
% Copyright 1996 - Ruben Caballero 
% New Mexico State University 
% 

% This function uses the mean and variance calculated by the function SPWSOBRE.M to 
% calculate the Average Symbol Variance. The result is a plot of the Average Symbol Variance 
% versus a for the SRRC filter as the spectrum shaping filter. 

clear; 

%* ********************* ******************************************** 

%** ********* * srrc Filters ****************************************** 

o/ o * *********** ****** ************************************************ 


alpha=[0. 1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1.0]; 
%SRRC = 1 (using delay=521, ascii file: samprcl.asc) 

meandeg(l,:) = [1.71 1 le+001 6.2204e+001 l.0710e+002 1.521 le+002 -1.6292e+002 
-1.1793e+002 -7.2849e+001 -2.7908e+001]; 

vardeg(l,:)= [3.2534e-003 3.0962e-003 3.7946e-003 3.8191e-003 3.0762e-003 3.4823e-003 
2.8293e-003 3.760 le-003]; 
avgvardegl = sum(vardeg(l,:))/8; 


%SRRC = 0.95 (using delay=521, ascii file: samprc6.asc) 

meandeg(9,:) = [1.7127e+001 6.2099e+001 1.0722e+002 1.5195e+002 -1.6292e+002 
-1.1 800e+002 -7.291 le+001 -2.7973e+001]; 

vardeg(9,:) = [9.2975e-003 1.2104e-002 8.6912e-003 1.2401e-002 1.5264e-002 1.0330e-002 
1.0367e-002 1.4794e-002]; 
avgvardeg9 = sum(vardeg(9,:))/8; 
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%SRRC = 0.90 (using delay=522, ascii file: samprc 1 9.asc) 

meandeg(19,:) = [1.7205e+001 6.2291e+00I 1.0734e+002 1.5186e+002 -1.6286e+002 
-l.I816e+002 -7.2892e+001 -2.805 le+001]; 

vardeg(]9,:) = [4.5507e-002 6.7076e-002 4.2139e-002 6.1049e-002 8.0428e-002 6 1238e-00? 

4.5872e-002 6.4998e-002]; 
avgvardeg!9 = sum(vardeg(19,:))/8; 




%SRRC — 0.85 (using delay=523, ascii file: samprc5.asc) 

meandeg(8,:)= [1.7292e+001 6.2435e+001 1.0743e+002 1.5174e+002 -1.6279e+002 
-1.1 834e+002 -7.2930e+00 1-2.8151 e+00 1 ]; 

vardeg(8,:)=[1.1523e-001 1.5848e-001 I.0640e-001 1.4940e-001 1.9194e-001 1 4397e-001 
1.1206e-001 1.4847e-001]; 
avgvardeg8 = sum(vardeg(8, :))/8; 

o/ Q ************** ********************** + * + ^ + + + + + 

o^*********** alpha=0.8 fs= 33 ***************************************** 

%* ******** * * * * * * * * * * * * * * * * * * * * * ♦ * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 


%SRRC - 0.8 (using delay=523, ascii file: samprc 1 8. asc) 


meandeg(18,:) = [1.7320e+001 6.2427e+001 1.0752e+002 1.5156e+002 -1.6276e+002 
-1.1 840e+002 -7.2974e+001 -2.824 le+001]; 

vardeg(18,:) = [1.5747e-001 1.9982e-001 1.4313e-001 1.9965e-001 2.6472e-001 1 8456e-001 
1.5906e-001 2.0486e-001]; 
avgvardeg!8 = sum(vardeg(18,:))/8; 


%** ********* *********************************** 

%*********** alpha=0.75 fs= 33 **************************************** 

%* ********** ******************************************************** 

%SRRC = 0.75 (using delay=524, ascii file: samprc4.asc) 

meandeg(7,:) = [1.7400e+00I 6.2585e+001 1.0763e+002 1.5153e+002 -1.6272e+002 
-1.185 le+002 -7.299 le+001 -2.8415e+001]; 

vardeg(7,:) = [2.522 le-001 3.0434e-001 2.4381e-001 2.9943e-001 4.1604e-001 2.9389e-001 
2.4814e-001 3 .0076e-00 1 ]; 
avgvardeg7 = sum(vardeg(7,:))/8; 
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%*** ******************** ********************************* *********** 
o^* ********** a iph a =o.70 fs= 33 **************************************** 
%* ***************************************************************** * 


%SRRC=0.7 (using delay= 524, ascii file: samprcl7.asc) 

meandeg( 1 7,:) = [ 1 ,7472e+00 1 6.2628e+00 1 1 ,0776e+002 1 .5 1 39e+002 - 1 ,6270e+002 
-1.1846e+002 -7.3051e+001 -2.8555e+001]; 

vardeg(17,:) = [3.0576e-001 3.5998e-001 3.0909e-001 3.5684e-001 5.0074e-001 3.5077e-001 
3.1454e-001 3.6826e-001]; 
avgvardegl7 = sum(vardeg(17,:))/8; 

%* ****************************************************************** 

%* ********** a iph a =0.65 fs= 33 **************************************** 
0 ^************************************************* * ***************** 


%SRRC = 0.65 (using delay=525, ascii file: samprc7.asc) 

meandeg(6,:) = [1.7547e+001 6.2769e+001 1.0778e+002 1.5138e+002 -1.6262e+002 
-1.1851e+002 -7.3048e+001 -2.8670e+001]; 

vardeg(6,:) = [4.0793e-001 4.6687e-001 4.1356e-001 4.4883e-001 6.4636e-001 4.5084e-001 
3.9264e-001 4.4991e-001]; 
avgvardeg6 = sum(vardeg(6,:))/8; 


%* ****************************************************************** 
o^*********** a iph a =o.60 fs= 33 **************************************** 
% * * ** * * **** **** * ** * **************** * **** *************** * ** * * * ****** * 


%SRRC=0.6 (using delays 525, ascii file: samprcl6.asc) 

meandeg(16,:) = [1.7568e+001 6.2800e+001 1.0775e+002 1.5123e+002 -1.6252e+002 
-1.1847e+002 -7.3137e+001 -2.8757e+001]; 

vardeg(16,:) = [4.81 62e-00 1 5.3695e-001 4.7133e-001 5.0882e-001 7.2892e-001 4.9668e-001 
4.4502e-001 5.1539e-001]; 
avgvardeg!6 = sum(vardeg(16,:))/8; 


%* ********** a ]ph a =o 55 fj= 33 **************************************** 


%SRRC = 0.55 (using delay=526, ascii file: samprc8.asc) 

meandeg(5,:) = [1.7492e+001 6.2998e+001 1.0763e+002 1.5124e+002 -1.6239e+002 
-1.1 850e+002 -7.3 1 1 4e+001 -2.8836e+00 1 ]; 

vardeg(5,:)= [5.6558e-001 6.0821e-001 5.2558e-001 5.8436e-001 8.3203e-001 5.5564e-001 
4.757 le-001 5.4309e-001]; 
avgvardeg5 = sum(vardeg(5,:))/8; 
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%SRRC — 0.5 (using delay — 529, ascii file: samprc2.asc) 

meandeg(4,:) = [1.7304e+001 6.3239e+001 1.0737e+002 1.5178e+002 -1.6229e+002 
-1.1858e+002 -7.277 le+001 -2.8588e+001]; 

vardeg(4,:) = [6.5188e-001 6.0204e-001 5.4793e-001 6.7353e-001 8.6404e-001 6 0872e-001 
4.6567e-001 4.82 12e-00 1 ]; 
avgvardeg4 = sum(vardeg(4,:))/8; 




%SRRC — 0.45 (using delay — 530, ascii file: samprc9.asc) 

meandeg(3,:) = [1.7206e+001 6.3302e+001 1.0730e+002 1.5201e+002 -1.6223e+002 
-1.1857e+002 -7.2514e+001 -2.8389e+001]; 

vardeg(3,:) = [7.25 1 1 e-00 1 5.9837e-001 5.5998e-001 6.7928e-001 8.3062e-001 6 4724e-001 
4.8553e-001 4.6707e-001]; 
avgvardeg3 = sum(vardeg(3,:))/8; 




%SRRC=0.4 (using delay= 529, ascii file: samprcl4.asc) 

meandeg(14,:)=[1.7222e+001 6.3314e+001 1.0752e+002 1.5 155e+002 -1.621 9e+002 
-1.1854e+002 -7.2540e+001 -2.858 le+001]; 

vardeg(14,:) = [7.4356e-001 6.0301e-001 5.5130e-001 6.2800e-001 8.3997e-001 6.3238e-001 
4.7616e-001 4.6947e-001]; 
avgvardegl4 = sum(vardeg(14,:))/8; 




%SRRC=0.35 (using delay= 530, ascii file: samprclO.asc) 

meandeg(2,:)=[1.7173e+001 6.3371e+001 1.0727e+002 1.5196e+002 -1.6212e+002 
-1.1 848e+002 -7.2149e+001 -2.8407e+001]; 

vardeg(2,:)= [8.1269e-001 5.5102e-001 5.7030e-001 5.9333e-001 7.6520e-001 63595e-001 
4.7928e-001 4.5532e-001]; 
avgvardeg2 = sum(vardeg(2,:))/8; 
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% :""‘~ 



%SRRC=0.30 (using delay= 529, ascii file: samprcl2.asc) 

meandeg(12 ) = [1.7161e+001 6.326 le+001 1.0737e+002 1 . 5159 e+ 002 -1.6217e+002 
meanaegi. ,-j l , 18 7 7e +o 02 -7 2458e+001 -2.8654e+001]; 

vardeg(12,:) = [7.2397e-001^5.3567e-001 5.4030e-001 5.1729e-001 7.4123e-001 5.5889e-001 
4.3710e-001 4.4035e-001]; 

avgvardegl2 = suni(vardeg(12,.))/8. 


% ************T , =n „ - „ 

- 

%SRRC =0.25 (using delay = 530, ascii file: samprclO.asc) 

meandeg(IO,)= [1.7306C+001 ^lOle+OOl 10709e4d02 1.5227e*»2 -1.6221**002 

® ’ i I813e-H)02 -7 2247e+001 -2.8346e-M)01]; 

vardeg(10,:)=[7.0668e-001 4.8405e-001 5.1696e-001 4.7003e-001 6.3451e-001 5.6028e-001 
4.2029e-001 3.7706e-001]; 
avgvardeglO = sum(vardeg(10,:))/8; 

..........*..*.**•*••*****•******** 

/o *** ,,****•****♦************************** 



%SRRC=0.2 (using delay= 529, ascii file: samprcl l.asc) 

meandeg(ll ■) = [1.7496e+001 6.2861e+001 1.0736e+002 1.5191e+002 -1.6257e+002 
meandegv ,.) 1 . 1796e +002 -7 2788e+001 -2.8484e+001]; 

vardeg(l 1,:) = [5.4128e-001 3.9687e-001 4.2458e-001 3.7001e-001 5.6274e-001 4.00524-001 
3.6781e-001 3.4022e-001]; 
avgvardegl 1 = sum(vardeg(l l,:))/8; 

/0 ., « 1S ; r _ -7! **************************************** 

* 

%SRRC=0.15 (using delay- 529. ascii file: sampiclS.asc) 

meandeg( 1 5 ) - [ 1 7699e+00 1 6.2560=+001 1.0740e+002 1.5213e+002 -1.6284e+002 

8 .1 I792e+002-72878e+001-2.8250e+001]; ... 

va,deg(15,:)- [4.0597e-001 3.2023e-001 3.5618e-001 2.8089e-001 4.1931e-001 2.9815e-001 
3.1036e-001 2.9837e-001]; 
avgvardegl 5 = sum(vardeg(15,:))/8, 




******************* 


******************* 


********** 
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**********************************************+******************* 
%*********** alpha=0.1 fs=33 ****************************************** 

o/ 0 ********* ********************** *********** ************************* 

%SRRC=0.1 (using delay= 529, ascii file: samprcl3.asc) 

meandeg(13,:) = [1.7752e+001 6.2291e+001 1.0743e+002 1.5202e+002 -1.629 le+002 
-1.1790e+002 -7.2870e+001 -2.81 13e+001]; 

vardeg(13,:) = [2.5572e-001 2.4821e-001 2.4306e-001 1.8168e-001 2.5692e-001 2 0489e-001 
1.9270e-001 2. 1 830e-00 1 ]; 
avgvardegl3 = sum(vardeg(13,:))/8; 


<y Q ************** pLQjg *************** 


% Average Variance for the Symbols 

figure(l); 

whitebg; 

AMPMavgvar = [avgvardegl3 avgvardegl5 avgvardegl 1 avgvardeglO avgvardegl2 avgvardeg2 
avgvardegM avgvardeg3 avgvardeg4 avgvardeg5 avgvardegl 6 avgvardeg6 
avgvardegl 7 avgvardeg7 avgvardegl 8 avgvardeg8 avgvardegl 9 avgvardeg9 
avgvardegl]; 

% interpolation using spline fit 

interaipha=0. 1 :0.000 1:1; 
interpol=spIine(alpha,AMPMavgvar,interalpha); 

pIot(alpha,AMPMavgvar/ok\interalpha,interpoI,'k’); 

title('SRRC Filter Average Symbol Variance vs Roll-of Factor (alpha) 1 ); 

xlabelfSRRC Rolloff factor (alpha)’); 

ylabel(' Average Variance (degrees)’); 

grid; 


201 


t 




